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In this paper, we investigate the testing problem that the spectral density matrices of several, 
not necessarily independent, stationary processes are equal. Based on an L 2 -type test statistic, 
we propose a new nonparametric approach, where the critical values of the tests are calculated 
with the help of randomization methods. We analyze asymptotic exactness and consistency of 
these randomization tests and show in simulation studies that the new procedures posses very 
good size and power characteristics. 
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1. Introduction and motivation 

Suppose, one observes two stretches of time series data Xi,..., and Yi,..., Y„ gen¬ 
erated by some unknown stationary time series processes. The question whether the 
models where these two stretches come from are identical is of considerable interest and 
has been investigated elaborately in statistical literature. For instance, by imposing cer¬ 
tain parametric assumptions on the data generation process (DGP), it becomes possible 
to compare suitable estimated parameters to judge how “near” the related models are. 
The assumption of underlying autoregressive processes of some hnite order, say, allows to 
compare estimated coefficients to quantify the nearness of both DGPs. But these models 
are usually not valid if the parametric assumption is not fulfilled. Therefore, nonparamet¬ 
ric methods are desired that can be justihed theoretically for a broader class of stochastic 
processes. 

If one is interested in second order properties, it seems obvious to focus on autoco¬ 
variances 7 (/i), autocorrelations p{h) or spectral densities /(w), where the latter choice 
has some advantages over the others, see, for example, Paparoditis [21]. Over the years. 
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several test statistics have been introduced that could be used for testing the null hy¬ 
pothesis 


{/jf(w) = /v(w) for all w e [-7t,7t]} (1.1) 

versus the alternative {fx^fv on a set of positive A\-measure}. These may be perio- 
dogram-based, see, for example, Coates and Diggle [6], Potscher and Reschenhofer [24], 
Diggle and Fisher [10], Caiado et al. [4, 5] and Preuss and Hildebrand! [25], or based on 
kernel spectral density estimators, see, for example, Maharaj [20], Eichler [11] or Dette 
and Paparoditis [9]. Jentsch and Pauly [16] have investigated the asymptotic properties 
of test statistics closely related to “periodogram-based distances” proposed by Caiado 
et al. [5]. They show that contrary to integrated periodograms these tests are usually 
not consistent which is caused by the use of non-consistent estimators for the spectral 
densities, the periodograms. In contrast to that, L 2 -type statistics based on the smoothed 
periodogram are shown to be consistent, see Paparoditis [21, 22], Eichler [11], Dette and 
Paparoditis [9] and Jentsch [15]. A possible test statistic for testing equality of spectral 
densities is therefore a suitable inflated version of 

[ (7x(w)-7y(a;))^dw, (1.2) 

J — n 


where fx and fy are nonparametric kernel estimators of the spectras. Central limit 
theorems for such statistics can be found in the above mentioned literature. Por instance, 
Paparoditis [21] has emphasized that statistics of this kind converge very slowly to a 
normal distribution, which results in a poor performance of the corresponding asymptotic 
test, particularly for small sample sizes. Therefore, other approaches as bootstrap or 
resampling techniques in general are needed to improve the small sample behavior of 
these tests. Paparoditis [21, 22] has used a (parametric) ARMA-bootstrap for goodness- 
of-fit testing and Dette and Paparoditis [9] have proposed a different approach based on 
the asymptotic Wishart-distribution of periodogram matrices. 

In this paper, we consider the more general testing problem of comparing the spectral 
density matrices of several, not necessarily independent, stationary processes. As moti¬ 
vated above, we will use an L 2 -type statistic of the form (1.2) which is adjusted for this 
multiple sample testing problem. To overcome the slow convergence speed of our test 
statistic, we propose a resampling method to compute critical values, the randomization 
technique, that exploits the well-known asymptotic independence of periodograms /„(a;) 
and /„(A) for different frequencies w ^ A with w, A G [0,7 t]. A nice feature of this ap¬ 
proach is that only one tuning parameter, the bandwidth of the involved kernel spectral 
densities, has to be selected. This can be done automatically by selection procedures 
based on cross validation as has been proposed by Beltrao and Bloomfield [2] and Robin¬ 
son [28]. In comparison to other mainstream resampling techniques as for example, the 
autoregressive sieve or block bootstrap techniques, it has the advantage that no other 
tuning parameter as order or block length has to be assessed in addition to the band¬ 
width. It is worth noting that these choices appear to be very crucial in applications and 
the performance of the corresponding tests usually reacts very sensitive on the choice 
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of these quantities. Nevertheless, these procedures are not applicable (at least) directly, 
because they do not mimic the desired null distribution under the alternative automat¬ 
ically. Moreover, the usage of randomization methods is very natural and has also been 
proposed by Diggle and Fisher [10] for the special (and included) testing problem of com¬ 
paring the spectral densities of two independent univariate stationary time series. This 
idea has been adopted by Maharaj [20] for comparing evolutionary spectra of univariate 
nonstationary time series. In both papers, the authors (only) account for the application 
of this procedure with the help of simulations and the following nice heuristic. Recall 
that a randomization test holds the prescribed level a exact for finite sample sizes if 
the data is invariant under the corresponding randomization group, see, for example, 
Lehmann and Romano [18]. Since In,x{^) and Inpf(uj) are asymptotically exchangeable 
in the case of independent time series, one would suggest that a randomization test, 
say , for (1.1) based on periodograms or functions of periodograms is at least asymp¬ 
totically exact, that is, E(ip’^) a holds as n —?► oo. However, this is only a conjecture 
and so far the randomization technique has not been analyzed in this context from a 
mathematical point of view in the literature. This is without question meaningful since 
it is offhand not clear in which situations it is applicable. We will close this gap in the 
next sections by investigating more general and in mathematical detail whether at all 
and, if true, under which possibly needed additional assumptions randomization-based 
testing procedures lead to success, that is, asymptotically exact and consistent tests. 
Amongst others it will turn out that this can be done effectively for testing the special 
null hypothesis (1.1). 

The paper is organized as follows. In the next Section 1.1, we introduce our notation, 
state the general testing problem and define the L 2 -type test statistic. Section 2 presents 
the main assumptions on the process and states asymptotic results for the test statistic 
which lead to an asymptotically exact and consistent benchmark test. The randomization 
procedure together with a corresponding discussion about its consistency is described in 
Section 3. Based on our theoretical results we compare the proposed tests in the special 
situation of testing for (1.1) in an extensive simulation study in Section 4, where the 
effects of dependence and bandwidth choice on the size and power behavior of both 
types of tests are investigated. All proofs are deferred to Section 6. 

1.1. Notation and formnlation of the problem 

We consider a d-dimensional zero mean stationary process (Aj, t £ Z) with d = pg, p, g £ N 
and g > 2. Under suitable assumptions, the process (Xf,t £ Z) posses a continuous (d x d) 
spectral density matrix f given by 



h— — co 


OJ £ [—7T,7T] 
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where r(/i) := E(Xf , h G Z, are the corresponding autocovariance matrices. To 


introduce our testing problem of interest, we write 



( /iiM 

/l2(w) •• 

■ /id(w)\ 



/2l(w) 

/22(w) 

/2d(w) 


f(w) = 






V/di(w) 

/d2(w) •• 

• fddipj)^ 



/Fn(w) 

Fi2(a;) 

■■ Fi,(w)\ 

= 

F2i(w) 

F22(w) 

F2q(w) 



VFgl(w) 

Fg2(w) • 

•• F,,(w)/ 


where Fmra(w) are {p x p) block matrices with Fmra(w) = F„m(w) for all m,n. Here 
and throughout the paper, all matrix-valued quantities are written as bold letters, all 
vector-valued quantities are underlined, A denotes the complex conjugation and the 
transposition of a complex matrix A. 

In this general setup, we want to test whether all {p x p) spectral density matrices Ffcfc 
of the q (sub-)processes {2Lt,kT^^ fc = 1, ■ • ■ ,9 with 


^t,k — (.^t,(k-i)p+i ,■■■, Xt^kp)^ 

are identical. Precisely, suppose we have observations ^ X.n hand and we want 

to test the null hypothesis 


Ho'. {Fii(a;) = F 22 (w) = •• • = F,,(a;) for all w G [-7t,7T]} (1.3) 

against the alternative 

Hi'. {3ki,k2 G{l,...,g},AcH with A\(A) > 0 such that Ffcjfcj(w) ^ Fk 2 k 2 {^) Vw G A}, 


where denotes the one-dimensional Lebesgue-measure on the Borel cr-algebra B. Ob¬ 
serve that this general framework includes (1.1). 

- T 

Consider now the periodogram matrix I(a;) := J(w) J(a;) based on JLi, ■ ■ ■ ,2Lni where 

1 ” 

^(‘^) •=(1-4) 

V27Tn 


is the corresponding d-variate discrete Fourier transform (DFT). For a better comparison 
with the spectral density matrix f, we write as above 


( /ii(a;) 

4i2(w) • 


d2i(w) 

422(w) 

hdi^) 

V4di(w) 

4d2(w) • 

■ Iddi^) J 
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/iiiM 

Il2(w) • 

• Il9(w)\ 

l2l(w) 

l22(w) 

l2g(w) 


1,2 (w) ■ 



where Imn(w) are (p x p) block matrices with Im„(a;) = for all m,n. Moreover, 

we define its pooled block diagonal periodogram matrix by 

1 9 

I(a;fc) := 

^ i=i 


With that we can introduce the kernel estimator 


^ 1 

f(a;) := - Kh{u} - u;k)liu}k), 

n ^' 

fe=-L(n-l)/2j 


W G [—7t, 7t] 


for the spectral density matrix f(aj), where [xj is the integer part of x G K, Wfc := 
27t^, fc = — ,..., are the Fourier frequencies, if is a nonnegative symmetric 

kernel function satisfying f if (x) dx = 27T, h is the bandwidth and Kh{-) := Note 

that the periodogram is understood as a 27T-periodically extended function on the real 
line throughout the paper. Moreover, we write as above 


^7ii(w) 

7i2(w) 

721 (w) 

722 (w) 

v.7di(w) 

fd2{uj) 

/Fii(w) 

Fi2(w) 

F2i(a;) 

F22(w) 

\F,i(w) 

Fg2(w) 


/id(w)^ 

hd{uj) 

fddiuj)/ 

Fi,M\ 

F2q(a;) 

Fgg(w)/ 


where Fmn(w) are uniformly consistent estimators of the {p x p) block matrices Fmn(w) 
for all m, n, that is, convergence sup^ ||Fmn(w) — Fm„(a;)|| —>■ 0 holds almost surely. Here 
and throughout the paper, || • || denotes the Frobenius norm, that is, for a matrix A = 
iaij)i<i,j<p G CP^P we set 


II Af :=tr(AA'^) = 

ij = l 

where tr(-) stands for the trace of a matrix, i.e. tr(A) := 
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For testing Hq, we now propose the L 2 -type test statistic 



(1.5) 


n 





where F(a;) := i the pooled spectral density estimator. 

For a better understanding of the asymptotic results in the following sections, it is 
interesting to note that 



holds, where Sjr = 1 ii j = r and 6jr = 0 otherwise. 

2. The unconditional test 

For deriving an asymptotically exact test for (1.3), we will in the sequel apply a result of 
Eichler [11] for a general class of L 2 -type statistics that includes the form (1.5). Therefore, 
we have to impose some assumptions. 

2.1. Assumptions 

First, we assume some strong mixing conditions for the DGP (compare Assumption 3.1 
in Eichler [11]). 

Assumption 2.1. (X^.t 6 Z) is a zero mean d-variate (strictly) stationary stochastic 
process defined on a probability space (fl, P). Furthermore, for any fc > 0, the kth order 
cumulants of {2Lt,t G'Z) satisfy the mixing conditions 






(ui,...,Ufc_i)l <oo 


( 2 . 1 ) 




for all j = 1,... ,k — 1 and oi,...,aj, = 1,...,d, where 


(ui,...,Mfc_i) =cum(A 






is the kth order joint cumulant of ■ ■ ■ ,Xuk_i,ak-n^o,ak (cf Brillinger [3] for the 

definition). 
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Note that the above assumption requires the existence of the moments of all orders. 
However, in the case that we presume that our process has a linear structure, it can be 
weakened. See Remark 3.2 below. 

Our second assumption on the kernel function K and the bandwidth h ensures the 
consistency of the kernel density estimators. It is similar to and implies Assumption 3.3 
of Eichler [11]. 


Assumption 2.2. 

(i) The kernel K is a bounded, symmetric, nonnegative and Lipschitz-continuous func¬ 
tion with compact support [—tt, tt] and 

/ Tt 

K{uj) dw = 271. 

-7T 

Furthermore, K(uj) has continuous Fourier transform k(u) such that 

J k‘^{u)du<oo and J k'^{u)du<oo. (2.2) 

(ii) The bandwidth h = h{n) is such that h^^^n —>■ 0 and h^n —^ oo as n —>■ oo. 


In comparison to the notation in Eichler [11], remark that a factor 27T is incorporated 
in K. Moreover, for a better lucidity we define the positive constants 

Ak ■= — J K^{v)dv and Bk'■= — J K(v)K(yz)dv'^ dz. (2.3) 

Furthermore, we like to point out that the optimal rate oi h (by means of an 

averaged mean integrated squared error type criterion) for estimating the spectral density 
is not covered by our assumptions. However, assumptions of such kind are commonly 
imposed in the literature in order to reduce the bias of the smoothed kernel estimator 
leading to a central limit theorem for r„, see Theorem 2.1 below. Compare for instance 
Taniguchi and Hondo [30], Taniguchi et al. [31], Taniguchi and Kakizawa [29], Eichler 
[11] or Dette and Paparoditis [9], who used similar (non-optimal) assumptions on the 
rate of the bandwidth. 


2.2. Asymptotic results for 

We are now ready to state a CLT for the test statistic T„ given in (1.5). 

Theorem 2.1 (Asymptotic null distribution of r„). Suppose that Assumptions 2.1 
and 2.2 are fulfilled. If Hq is true, it holds 
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as n —?► oo, where 


1^0 = Ak 





(2.4) 


and 


L 


■n 


Bk 


— rc \ Q 


^ X! (“ 1 + 9^jih ) (“ 1 + ) 


jl,j2 J3 J4 = l 



(2.5) 


where Fjj(oj) = Fn (w) holds for all lj and j = 1,..., g. Here and throughout the paper, 

'D 

—> denotes convergence in distribution. 

Remark that we have chosen the above representation of and Tq for notational 
convenience. It is of course possible to split the sums above and to rewrite the expressions 
by using = Fn. 

It is interesting to note that we do not need the assumption of positive definite spec¬ 
tral density block matrices since we do not work with ratio-type test statistics, see, for 
example, Eichler [II] for p = I and g = 2 or Dette and Paparoditis [9] for the case p = I. 

We give one important example which is possibly of most relevant interest, where we 
compare two real-valued spectral densities. 

Example 2.1 (The case p = I and q = 2). For p = 1 and g = 2, the quantities po and 
Tq defined in Theorem 2.1 become 



rc 


and 



where Cjk{u}) = /- ^ ifjji^)fkk (^) > 0) is the squared coherence between the two 

components j and k of {2(_t,t G Z) (cf. Hannan [12], page 43). 




Since po and Tq are in general unknown we have to estimate them before we can apply 
the above results for testing (1.3). Following Eichler [11], Remark 3.7, we can estimate 
both quantities by substituting Fjk in their expressions by the corresponding consistent 
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estimators Fjk- Rewriting the right-hand side of the equations (2.4) and (2.5) under Hq 
by means of the equality Fjj = ^ J2i=i leads to 

^■=Ak [ [(g-l)|tr(F(a;))|^ - i ^ |tr(Fj,j,(w))|^^ dw (2.6) 

^n,j2=i ^ 

3i¥^h 

as an adequate estimator for fio. Similarly, we construct a consistent estimator for Tq . 

Note that it is in applications more convenient to use a discretized version of /r, where 
the last integral is approximated by its Riemann sum. 

Now we are ready to define the unconditional test 

Tn l(tii_a,oo) ((^n h ^ 

where rti_a denotes the (1 — a)-quantile of the standard normal distribution. Its prop¬ 
erties are summarized in the following theorem. 

Theorem 2.2. Suppose that the Assumptions 2.1 and 2.2 hold. Then the test tpn is not 
only of asymptotic level a, that is, E{ipn) —> a holds if Hq is true, but also consistent for 
testing Hq versus Hi, that is, E((pn) —> 1 under Hi as n —^ oo. 

2.3. Power under local alternatives 

When studying the behavior of the test ipn under local alternatives we have to consider 
observations Xf.... , X” that come from sequences of d-dimensional zero-mean processes 
{Xff,t G Z) with spectral density matrices given by 

f”(w) = f(w)-f a„g(w). (2.7) 

Following Eichler [11] we suppose that f is nonnegative definite, Hermitian and sat¬ 
isfies Hq and that —>■ 0 as n —>■ oo. Moreover, f and g are assumed to be twice 
continuously differentiable. In addition, for every u G N the sequence of processes 
{Xff ,t G Z) satisfies the strong mixing condition (2.1) uniformly in n, where we replace 
Cai,....afc(wi,---,Ufe_i) by cum(X)^^ (the fcth order joint cumu- 

lant of Under these assumptions it can be deduced from 

Theorem 5.4 in Eichler [ll] that the test can detect local alternatives up to an order 
of To be concrete, we have the following result. 

Theorem 2.3. Under the conditions stated above the asymptotic power of the test ipn 
under local alternatives (2.7) with an = is given by liningoo E{ipn) = 1 ~ 

4>(ui_q, — u/tq), where $ denotes the c.d.f. of the standard normal distribution and the 
detection shift v is defined as 

r — T 

v= j vec(g(a;)) Fvec(g(a;))dtu. 


— 7T 
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Here the {d^ x d^) matrix F = (rr,s)r,s has the entries 

for \i-k\ = \j -l\=0 and {i,j) G S, 

for \i — k\ = \j — l\€ pN and {i,j) G S, 
otherwise, 

with S = lJc=i{(^’'®)’ 1 + (c — l)p <r,s< cp}. 

Note, that the same detection order has also been found by Paparoditis [21]. 


- i-\-{j—l)d,k-\-{£—l)d 


= < 


g-i 

q 
1 


lo 


3. The randomization tests 

Although the asymptotic performance of the unconditional test proposed above seems 
to be satisfactory, it is well known that the convergence speed of L 2 -type statistics like 
(1.5) is rather slow (see Hardle and Mammen [13] as well as Paparoditis ([21, 22]) for 
more details). For small or moderate sample sizes, we therefore expect that the above 
Gaussian approximation will not perform pretty well in applications. Our simulations in 
Section 4 support this presumption. To overcome this hitch, we propose in the following 
a suitable randomization technique for approximating the finite sample distribution of 
our test statistic. This method leads to so called randomization tests that use more 
appropriate, data-dependent critical values. Let tt^ : (f2, JF, IP) —5,, fe = 0,..., , with 

7rfe = (7rfe(l),...,7rfc(q))^ (3.1) 

be a sequence of independent and uniformly distributed random variables on the symmet¬ 
ric group Sg (the set of all permutations of (1,..., q)) dehned on some further probability 
space (n, J", P). In what follows, we assume that {Trk)k and are independent (defined 
as random variables on the joint probability space (fl x id, if G P)). Conditional 

on the observations 2Li j ■ • ■ j 2Ln fhe corresponding randomization test statistic is given 

by 


rix 9 


Tl=nh^l^ \ Y. 


Ln/2J 


^ ^ U;/2)(l7rfc (r),7rfc(r) (^A:) I(^A:)) 


fc=-L(n-l)/2j 


dw, 


where we set 7r_fc := and TTk+sn ■ = for s G Z to maintain the symmetry properties 

of a spectral density matrix. Note that similar to the unconditional case, it holds 

1 

l7rfc(r),7rfc(r) (Wfc) — I(Wfc) = —— ^^(1 ~ (r-) (Iji ip>k)- 

^ i=i 
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3.1. Asymptotic results for T* 


In the following, we will analyze in which situations our randomization method leads to 
asymptotically valid tests, that is, whether some suitable centered T* converges to the 
same distribution as Tn — Therefore, we have to exploit the limiting behavior 

of the randomization statistic T*. It will turn out that we have to center T* at 


rii i (-1+9%..) 


h,h=i 


(3.2) 




to gain asymptotic normality. 


Theorem 3.1 (Asymptotic distribution of T*). Suppose that the mixing condi¬ 
tion (2.1) holds for all k <32. Under the Assumption 2.2, we have (conditioned on 
X ,,.... X„ ) convergence in distribution 


T:-^AAr(0,r*2) 


(3.3) 


in probability as n —>■ oo, where fi* is as in (3.2) and r*^ is defined by 

9 


r*' := Bk 


f f 1 + + 1 (1 '^ili3)(l %j4)l 

.1..203J4 = % ^ ^ 


X {|tr(Fjjjg (u;)Fj2j4 (w) )| 


(3.4) 


+ tr(Fjiji (a;)Fj2j2 (^) ) f^’(l'i3i3 (^)%4i4 )}'lW' 


Note, that the conditions of the above theorem are weaker than the assumptions of 
Theorem 2.1. Moreover, remark that the centering term (2* defined in (3.2) fulfills (2* = 
p* + op^y/h), where 


P* '■= Ak 



X {|tr(Fj4j2(w))| +tiiFj^j^{uj)Fj^j^{uj) )} dw. 


(3.5) 


Hence, to compare the above results with Theorem 2.1 we have to analyze p* and r*^ 
under Hq which leads to p^ and Tq^ in the following remark. 
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Remark 3.1. Under Hq the constants fi* and of Theorem 3.1, reduce to 


Mo 


= Ak [ (- (-l + 9'^nt2){|tr(Fjij2(w))|^}] da; = 


Mo 


(3.6) 


and 


To*" = Bk 


f ri^ f 1 +'7'^iii3^i2i4 + 1 (1 ‘^ 1214 )') 

nj2j3j4=l'^ ^ ^ 


(3.7) 


X |tr(Fj,j3(w)Fj,j,(w) )| dw, 


where Fjj(aj) = Fii(a;) holds for all a; and j = 1,..., g. 


It is interesting to note that the unconditional and conditional centering parts gio in 
(2.4) and \x*q in (3.6) coincide. Moreover, T„ — as well as T* — posses a 

Gaussian limit distribution under the null. However, in order to construct an asymptot¬ 
ically exact randomization test based on T* — we have to be sure that the limit 

variances Tq in (2.5) and Tq^ in (3.7) under Hq are equal as well. This will be discussed 
in more detail in the next subsection. Here we just state an example (corresponding to 
Example 2.1) where this requirement is fulfilled. 


Example 3.1 (The case p=l and q = 2). For p = 1 and q = 2 the quantities p* and 
defined in Theorem 3.1 become 


P*=Ak j i{(/ii(u;)-/22(w))^-h/n(u;)-h/|2(w)-2|/i2(a;)|^}dw 


and 


= Bk 


^{(/ll(w) - f22{pj))‘^ 


(/ii(‘^) +/22('^) “ 21 / 12 ( 0 ;)I ) }da;. 


Under Hq, we have 


/ TT 

/ii(w)(l - Ci2{uj)) do; = pq 

-n 


and 


Tq‘^=Bk [ Ci2{uj)f duj = T^. 

J — 71 
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Based on the conditional CLTs above we can now define different randomization tests. 
The first natural approach is to use the test </3jj_cent •= ~ where 

c* cent(Q^) fhe data-dependent (1 — Q;)-quantile of the conditional distribution of T* — 
h~^/^fi* given the data. As typical for resampling methods, note that we still use the same 
test statistic as for the unconditional case and only apply the randomization statistic to 
calculate the critical value. By Theorems 2.1 and 3.1, this test will be asymptotically 
exact, i.e. —>■ a holds under Hq, if the asymptotic variances Tq and Tq^ of 

the test statistic T„ — h~^^'^fi and its randomization version posses the 

same limit under the null. Although T* — does in general not mimic the null 

distribution under the alternative, the following corollary shows that be 

asymptotically consistent in these situations. Before we state these properties, we like to 
introduce a computational less demanding version of the above test. Therefore note that 
Remark 3.1 implies that the difference — fi) converges to zero in probability 

under Hq. Hence, it may be convenient to use the test 

V’ri l(c* (q:),oo) (Ti) (3-8) 

without centering part, where c* (a) is the data-dependent (1 — a)-quantile of the condi¬ 
tional distribution of T*. For completeness, we shortly explain the numerical algorithm 
for the implementation of 1 , 5 *. The algorithm for 1 ^* is analogue. 

Step 1: Compute the test statistic T„ as given in (1.5) based on data A^,..., A„. 

Step 2: Generate independent random permutations ttq, ... ,7r[„/2j as in (3.1) that are 
uniformly distributed on the symmetric group Sq. 

Step 3: Calculate the randomization statistic T* (given 2Li^ ■ ■ ■ ■ ). 

Step 4: Repeat the Steps 2 and 3 H-times, where B is large, which leads to 

7 ■ • ■ 5 n 

Step 5: Reject the null hypothesis Hq if B~^ Y^b=i > a. 

In the sequel, we analyze the asymptotic properties of both tests and tp* and 
compare it with the unconditional benchmark test 

Corollary 3.1 (Exactness and consistency of and (/Jjljcent)- Suppose the as¬ 
sumptions of Theorem 2.1 are satisfied. 

(a) If Hq is true, the following assertions are equivalent: 

(i) The randomization test v?* is asymptotically exact and equivalent to ipn, that is, 

if(|(/j„ — yj* I) —>• 0 asn^oo. (3.9) 


0= ^ ((g-l)3-l)|tr(F,,,3(w)F,,,3(w)^)|^ 

flj 3 = l 


(ii) It holds 
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-2 ^ + (3.10) 

il J3 J4 = l 

all^ 

9 _ j. 

+ Y1 (9-2)|tr(Fj,j3(a;)Fj,j^(a;) )|^ 

3i,n,33 j4=i 
ii ¥=h J2 5^i4 ,ii ¥^32,33 ¥ 34 , 

Moreover, if (3.10) holds, will also be asymptotically consistent, that is, E(lp)() —>• 1 
under Hi asn^oo. 

(b) The above statement (a) also holds true for instead of ip^^. 

Note that the asymptotic equivalence in (a) implies that both tests posses the same 
power for contiguous alternatives. Moreover, remark that the above stated consistency of 
the tests is caused by the fact that r„ — h~^/‘^'jl converges to +00 in probability for fixed 
alternatives, see the proof section for more details, and the non-degenerated limit law 
of its randomization version. In the following, we will give some necessary and sufficient 
conditions for (ii) above. 

Corollary 3.2 (Necessary and sufficient conditions for exactness of and 

Tncent)- Suppose that the assumptions of Theorem 2.1 hold. 

(i) In the case q = 2, that is, we are testing for the equality of two (p x p) spectral 
density matrices, condition (ii) of Corollary 3.1 is fulfilled for all p G N. 

(ii) For any q,p £ N, q>2, a sufficient condition for condition (ii) of Corollary 3.1 
is 


Fy(a;) = Fi2(w) (3.11) 

for all uj £ [—7T,tt] and all i,j £ {1 ,... ,q} with i ^ j. This means that all (p x p) 
block matrices on all secondary diagonals have to be equal and Hermitian every¬ 
where. 

(iii) For p = l and q>2, a sufficient condition weaker than (ii) above is 

\f.,{u3)f = \fl2{u3)f 

for all UJ £ [—7T,7t] and all i,j £ {1, ... ,q} with i j. 

(iv) For p = 1 and q = 3, the condition 

|/4,(u;)|^ = |/i2(cu)r 

for all UJ £ [— tt, tt] and all i,j £ {1,2,3} with if=j is not only sufficient, but also 
necessary for condition (ii) of Corollary 3.1. 

Note that the randomization tests tp* and tp* (.gjjj, are asymptotically exact particu¬ 
larly in the case of uncorrelated p-variate (sub-)time series X_i,.. - due to (ii) above. 
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Intuitively, this makes sense because permuting the block diagonal matrices distorts the 
correlation structure between these time series, if there is any. This explains also why both 
randomization tests are asymptotically exact in the more general case of equal covariance 
structure between the time series X.ij ■ ■ ■ as shown in (ii) and (hi) of Corollary 3.2 
above, which is often denoted by equicorrelatedness. 

Observe that due to the sophisticated condition in Corollary 3.1, a result for general 
q>3 equivalent to (iv) in Corollary 3.2 does not seem to be achievable. 

Nevertheless, there is an obvious possibility for constructing another randomization 
test that works even if the presumption (3.10) of Corollary 3.1 is not fulfilled. We just 
have to impose an estimator for the conditional limit variance of T* as introduced 
in Theorem 3.1. An appropriate candidate is given by 


:= Bk 


l-n r 


9 

E 

jl,j2 J3 J4 = l 


1 + qSjij^Sj^j^ 


9-1 


(1 '^iii3)(l 


X {tr(Fjjjj (a;)Fj2j2 (uj) ) tr(Fj^j^(u’)Fj^j^(u)) ) (3.12) 

and the corresponding randomization test is stud l(c* * u(a),oo)((Tn — /t), 

where c* stud('^) denotes the data-dependent (1 — a)-quantile of the conditional distribu¬ 
tion of (T* — given the data. Let us shortly state its properties. 


Corollary 3.3. Suppose that the assumptions of Theorem 2.1 hold. Then the test (/?* 
is not only of asymptotic level a under Hq but also consistent for testing Hq versus Hi, 
that is, E{(p^ stud) —^ 1 under Hi as oo. 

Finally, we like to note that some of the assumptions can be weakened for linear 
processes. In particular, by combining Theorem 3.1 above with results of Dette and 
Paparoditis [9] we gain the following remark. 

Remark 3.2 (Exactness and consistency of the tests for linear processes). If 
the process posses a linear structure, that is, X_f = ^jSu-j for ^ d-dimensional 

i.i.d. white noise {e^,t G Z) and a sequence of {d x d)-matrices (^j)j = (('0i(r,s))r,s)j, 
Condition (2.2) is not needed and the mixing Assumption 2.1 in the above Corollaries 
3.1-3.3 can be substituted by the summable condition s)| < oo and the 

moment condition Adletp^) < oo. 


For the two sample testing problem as described in the introductory part our ran¬ 
domization test has a nice reading as a symmetry test. This is part of the following 
remark. 


Remark 3.3 (The case q = 2: Interpretation as a conditional symmetry test). 
Remark that we can rewrite (in distributional equality) the summands of our random- 
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ization statistic for q = 2 as 

~ "D 1 

l7rfc(r).7rfc(r)(Wfe) “ I(Wfc) = -efc(Iii(u;fc) “ l22(wfe)), 

where {ek)k>o are i.i.d. signs, that is, independent and on {+1, —1} uniformly distributed 
r.v.s, and we set e_fc := Ck as well as Ck+sn ■= Cfc for s G Z. Hence, the above defined 
randomization tests can be interpreted as some general kind of conditional symmetry 
tests in this sitnation. 

Remark 3-4 (Choice and influence of the bandwidth). As already mentioned in 
the Introduction, our randomization approach has the nice advantage that the bandwidth 
is the only tuning parameter that has to be assessed. This feature becomes even better 
as we will see in our extensive simulation study, see Section 4 below, where our approach 
does not react very sensitive to variations of the bandwidth h. However, it is of course 
desirable to have a detached principle for selecting the bandwidth. Therefore, we use 
a data driven cross validation method for choosing h for the kernel spectral density 
estimation. This method is due to Beltrao and Bloomfield [2] and Robinson [28] and has 
also been applied by Paparoditis [21-23]. 

Finally, we also like to discuss the local power of all randomization tests. 

Corollary 3.4. Suppose that the assumptions of Theorem 2.3 hold. 

(i) The test stud same local power as the asymptotic test ipn, i.e. we 

have £'((/?* —>• 1 — — v/tq) under local alternatives (2.7) with a„ = 

h~^l'^n~^l'^, where v is given in Theorem 2.3. 

(ii) In addition, suppose that f in (2.7) satisfies (3.10). Then and 
posses the same power under local alternatives (2.7) with an = h~^l^n~^l'^. 


4. Simulation studies 

In this section, we illustrate the performance of the randomization procedure in compar¬ 
ison with the asymptotic (unconditional) benchmark test as described in the previous 
sections. 

For better lucidity, we thereby only analyze the finite sample behavior of the com¬ 
putational least-demanding randomization test (/?* as proposed in (3.8). The other two 
randomization procedures from Section 3 behave similar to or slightly worse than </?* . For 
more details, we refer the reader to our supplementary material (cf. Jentsch and Pauly 

[17]). 

4.1. The setup 

Suppose we observe bivariate time series data (A^ = {Xtp,Xt^ 2 )'^,t = 1,... ,n) and we 
want to test the null hypothesis Hq of equality of both corresponding one-dimensional 


Testing equality of spectral densities 


17 


spectral densities /i(w) and / 2 (w). In the setup of Section 1, this means q = 2, p=l and 
fj{uj) = Fjj{uj), j = 1,2 and we test 

Ho- {/i(w) = / 2 (w) for all w £ [-7t,7t]} 


against 


Hi: {3A C S([—7t,7t]) with A\(A) > 0 : fi{uj) ^ fiiyi) for all w £ A}. 

In the following, we consider data from several well-established time series models. In 
particular, our analysis includes Gaussian and non-Gaussian linear time series as well 
as non-linear time series models. The linear models under consideration cover moving 
average (MA) models and autoregressive (AR) models with innovations following Gaus¬ 
sian, logistic and double-exponential distributions, respectively. GARGH models, thresh¬ 
old AR (TAR) and random coefficient autoregressive (RGA) models are investigated to 
cover important classes of non-linear time series. 

Although GARCH processes are known to have power law tails (see Basrak et al. [1], 
Section 4) and, consequently, only moments up to some finite order exists, we include 
GARCH models to investigate the general performance of the randomization approach 
for processes that go beyond our Assumption 2.1. For the same reasons, we consider also 
RCA models in our simulation study. 

The performance of the randomization test in comparison to the unconditional 
benchmark test ipn is investigated under the null and under the alternative. 

For all models under consideration, we have generated T = 400 time series. For evalua¬ 
tion of the test statistic, the bandwidth has been chosen by cross validation as proposed 
in Remark 3.4 and is denoted by hcv- Further, we use the Bartlett-Priestley kernel (see 
Priestley [27], page 448) for which the constants in (2.3) become Ak = | and Bk = 
for this particular kernel function. For each time series, the test ipn has been executed 
with critical values from normal approximation as discussed in Section 2 and the ran¬ 
domization test (/?* as discussed in Section 3, where B = 300 randomization replications 
have been used. 


4.2. Analysis of the size 

To investigate the behavior of the tests under the null, we consider realizations from 
vector autoregressive models (VAR), vector moving-average models (VMA) to cover lin¬ 
ear time series and from GARCH, TAR and RCA models to cover non-linear cases. We 
consider data from the bivariate vector AR(1) model 

AR(1): A,=AA,_i + e*, t€Z, (4.1) 

where ~ (0, Si) is an independent and identically distributed (i.i.d.) bivariate white 
noise with covariance matrix Si, A is chosen from 
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and Si = Id is the unit matrix. The VAR model corresponding to coefficient matrix 
Ai is denoted by ARi. All models have been investigated for i.i.d. Gaussian innovations 
(cjjt G Z), whereas the most critical model AR^ has also been analyzed for logistic (with 
c.d.f. F{x) = (1 + exp(—x))“^) and double-exponential distributions (with p.d.f. f{x) = 
exp(—|a;|)/2) of the i.i.d. innovations G Z), respectively. Observe that due to the 
diagonal shape of all involved matrices Si and A^, we are dealing with two independent 
univariate time series here. Furthermore, we consider data from the bivariate vector 
MA(1) model 


MA(1); A(Be(_i-I-gj, t G Z, 


(4.2) 


where ~ (0, S 2 ) is an i.i.d. bivariate white noise with covariance matrix S 2 , B is chosen 
from 


Bi 


fOA 0.5\ 

1 ^ 0.5 0 . 1 ) ’ 


and 


B 2 


/0.5 0.5\ 
1^0.5 0.5 ) ’ 





(4.3) 


The VMA model corresponding to B^ is denoted by MAi. Again all models have been 
investigated for i.i.d. Gaussian innovations (e^jt G Z), whereas the model MA 3 has also 
been analyzed for logistic and double-exponential distributions of the innovations, re¬ 
spectively, as discussed above, but with covariance matrix S 2 . In this setting, we are 
dealing with two dependent time series whose marginal spectral densities are equal due 
to the symmetric shape of S 2 and B^. Note that the application of the randomization 
technique to the case of two dependent time series is justified by Corollary 3.2(i). Also, 
we investigate three different non-linear time series models. First, we consider bivariate 
data (X^ = (At_i, At_ 2 )^,i = 1 ,... ,n) from two independent, but identically distributed 
univariate GARCH(1,1) processes {Xt^i, t G Z}, i = 1,2, with 

GARCH(1,1): Xt^i = = w-|-^-|-j, f G Z, (4.4) 

where lo = 0.01, a = 0.1 and the coefficient b is chosen from 


bi = 0 . 2 , 62 = 0.3, 63 = 0.4. 

The corresponding models are denoted by GARCHi, i = 1,2,3. Further, two (centered) 
independent, but identically distributed univariate TAR(l) models {Xt^i-, t G Z}, z = 1,2, 
that follow the model equation 


( a(l)Xt-i^i + et^i, At_i^i<0, 

\^a(2)Xt-i^i + et^i, > 0 , 


TAR(l): Xt,i = 


f G Z 


(4.5) 
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with coefficients a= (a(l),a(2))' chosen from 

«! = (-0.2,0.!)^, 02 = (-0.3,0.2)^, 03 = (-0.4,0.3)^ 

are studied and denoted by TARi, i = 1,2,3. Also two independent, but identically dis¬ 
tributed univariate RCA(l) models {Xt^i,t € Z}, i = 1,2, that follow the model equation 

RCA(l): X-t^i = atXt-i + t € Z, (4.6) 

where at ^ A/'(0,cr^) are i.i.d. centered normally distributed random variables with stan¬ 
dard deviation a chosen from 


tTi=0.1, 02=0.2, 03=0.3 

are considered and denoted by RCAi, i = l,2, 3. For all non-linear models above, we have 
used (independent) standard normal i.i.d. white noise processes {et^i}. 

For nominal sizes a € {1%,5%, 10%} and sample sizes nG {50,100,200}, the corre¬ 
sponding results for all combinations are displayed in Tables 1-3. To check how sensitive 
the tests react on the bandwidth choice, we report the simulation results for bandwidths 
c- /lev and cG {0.5,1,1.5} to cover under-smoothing and over-smoothing with respect 
to the bandwidth hcv chosen via cross validation. 

To illustrate the slow convergence of the actual size of the unconditional test (pn for 
the different models under the null and to emphasize the need of resampling techniques 
to resolve this issue, we report its performance also for larger sample sizes in Table 4. 

4.3. Analysis of the power 

As can be seen in Tables 1-3, the asymptotic test tends to overreject the null hypothesis 
systematically in most situations and should not be applied, at least for the simulated 
sample sizes. In particular, for the autoregressive models where the actually achieved size 
is far too large compared to nominal size (see Tables 1-2), the unconditional benchmark 
test ifn cannot be judged. Moreover, even for larger sample sizes up to n = 2000 it does 
not keep the prescribed level satisfactorily, see Table 4. Thereforee, we present here only 
the power behavior of the randomization test (^*. However, we present a small power 
simulation study of ipn as well as all other randomization tests in the supplementary 
material (cf. Jentsch and Pauly [17]). It can be seen that there is actually no big difference 
in the power behavior (measured as achieved power) between the asymptotic test ipn and 
all randomization tests. 

To illustrate the behavior of under the alternative, that is for inequality of both 
spectral densities, we consider several models belonging to the same model classes that 
have already been considered above under the null. First, we consider realizations from 
the autoregressive model in (4.1). Here, A is chosen from 



to 
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Table 1. Actual size of and for nominal size a G {1%, 5%, 10%}, sample size n € {50,100,200}, bandwidth h = c - hcv for 
c€ {0.5,1,1.5}, autoregressive models AR1-AR3 and moving-average models MA1-MA3, all with Gaussian innovations 




ARi 






AR 2 






ARi 







n: 

50 


100 


200 


50 


100 


200 


50 


100 


200 






















a 

c 

^71 


’•pn 

Vn 

pn 

pn 

pn 

(Pn 

pn 

Vn 

pn 

pn 

pn 

(Pn 

pn 

(Pn 

pn 

pn 

1 

0.5 

5.5 

1.5 

8.5 

2.0 

6.5 

1.0 

10.0 

1.5 

7.5 

1.3 

7.3 

1.8 

17.8 

2.8 

13.0 

1.3 

11.3 

1.0 


1 

11.3 

1.3 

10.5 

1.0 

11.3 

0.8 

13.8 

1.5 

11.3 

2.3 

12.0 

1.3 

27.0 

2.0 

20.0 

1.5 

21.8 

1.8 


1.5 

9.8 

1.5 

17.0 

1.8 

10.5 

1.8 

16.8 

0.5 

13.3 

1.0 

11.5 

1.5 

36.3 

2.0 

26.5 

1.3 

24.0 

1.0 

5 

0.5 

10.5 

5.0 

14.5 

5.8 

12.3 

4.8 

14.3 

6.0 

14.8 

4.5 

12.3 

6.8 

27.3 

7.5 

20.8 

5.8 

17.5 

4.8 


1 

17.0 

8.3 

17.3 

6.3 

17.3 

4.5 

20.5 

5.5 

17.8 

6.3 

18.5 

7.0 

33.5 

6.5 

27.3 

5.3 

27.3 

6.3 


1.5 

18.0 

5.3 

21.8 

8.3 

17.5 

5.3 

23.5 

4.8 

18.5 

4.8 

17.5 

4.8 

42.5 

7.3 

33.8 

4.3 

30.8 

5.8 

10 

0.5 

15.0 

8.3 

20.3 

11.8 

18.3 

10.8 

20.0 

11.0 

19.0 

10.0 

16.3 

11.0 

33.0 

12.8 

28.8 

11.5 

22.5 

10.5 


1 

19.8 

12.3 

22.0 

12.0 

23.8 

11.3 

24.8 

11.5 

23.5 

11.5 

21.3 

12.8 

38.0 

13.3 

33.5 

9.8 

32.8 

11.5 


1.5 

22.3 

8.0 

27.3 

13.3 

21.8 

11.8 

27.3 

11.3 

24.3 

9.0 

23.8 

9.3 

48.0 

13.8 

41.0 

8.5 

34.3 

10.0 



MAi 






MA 2 






MAz 







n: 

50 


100 


200 


50 


100 


200 


50 


100 


200 


a 

c 



pn 


pn 

pn 

pn 

iPn 

pn 


pn 

pn 

pn 

(Pn 

pn 

(Pn 

pn 

pn 

1 

0.5 

4.0 

2.0 

4.0 

1.5 

4.8 

1.8 

6.0 

1.8 

5.8 

0.5 

4.0 

1.3 

6.5 

1.3 

5.3 

1.8 

5.3 

1.3 


1 

5.0 

1.5 

4.3 

2.0 

5.5 

1.5 

7.5 

1.8 

4.3 

0.5 

5.5 

1.8 

8.3 

2.0 

6.5 

1.0 

5.8 

1.8 


1.5 

3.3 

0.5 

4.0 

1.5 

3.3 

0.5 

9.5 

1.3 

6.8 

1.3 

7.0 

1.0 

11.5 

1.0 

8.5 

2.0 

11.0 

2.5 

5 

0.5 

7.8 

5.8 

8.5 

5.0 

8.5 

6.8 

11.0 

5.5 

10.0 

5.3 

11.3 

5.3 

11.0 

5.8 

10.5 

6.3 

8.8 

5.5 


1 

8.8 

5.0 

9.5 

4.8 

10.0 

6.5 

13.5 

7.3 

8.0 

4.8 

11.3 

6.0 

13.5 

7.3 

12.0 

6.5 

11.8 

5.5 


1.5 

8.3 

4.8 

7.5 

4.3 

9.5 

6.3 

14.3 

6.5 

11.0 

5.3 

11.5 

5.3 

17.8 

5.0 

13.5 

5.5 

17.0 

7.5 

10 

0.5 

13.0 

8.3 

13.5 

10.3 

12.3 

11.5 

13.8 

9.5 

16.0 

10.3 

15.8 

10.3 

16.0 

9.8 

16.5 

10.8 

12.8 

11.3 


1 

12.8 

11.3 

14.0 

11.5 

12.0 

11.5 

18.5 

12.8 

13.0 

11.0 

16.8 

12.3 

17.8 

11.0 

15.8 

11.5 

16.3 

9.5 


1.5 

11.5 

10.5 

12.8 

10.3 

13.3 

13.0 

18.0 

13.5 

14.5 

11.0 

15.3 

10.8 

23.0 

8.8 

18.0 

9.3 

20.8 

12.8 
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Table 2. Actual size of ipn and (p* for nominal size a € {1%,5%,10%}, sample size n £ 
{50,100,200}, bandwidth h = c-hcv for c € {0.5,1,1.5}, autoregressive model AR 3 and moving- 
average model MA 3 , both with logistic and double-exponential distribution of the innovations 


{et} 

a 

n: 

c 

AR 3 






MA 3 






50 


100 


200 


50 


100 


200 



Tn 


Tn 


Tn 

^71 

Tn 


Tn 

9^ri 


Logistic 

1 

0.5 

16.3 

3.3 

15.8 

2.5 

14.3 

2.0 

6.8 

2.8 

6.3 

1.5 

8.0 

2.0 



1 

33.8 

4.0 

23.3 

0.8 

21.0 

1.8 

9.8 

1.8 

9.0 

2.0 

10.0 

2.8 



1.5 

39.3 

3.3 

36.3 

2.0 

29.8 

2.5 

16.5 

4.8 

14.5 

1.8 

10.5 

2.5 


5 

0.5 

26.8 

6.5 

23.8 

7.3 

21.0 

7.0 

12.8 

9.3 

12.0 

6.8 

13.5 

8.3 



1 

40.8 

8.0 

30.5 

7.3 

27.5 

4.8 

17.0 

8.5 

16.5 

8.3 

17.0 

9.3 



1.5 

46.3 

7.5 

45.3 

7.5 

34.8 

8.8 

24.0 

11.0 

20.0 

7.8 

17.5 

6.5 


10 

0.5 

33.3 

13.5 

27.8 

11.8 

24.5 

12.8 

17.5 

13.3 

16.5 

13.5 

17.0 

15.3 



1 

45.5 

12.3 

35.0 

11.5 

33.0 

10.3 

19.8 

13.5 

21.5 

14.0 

21.8 

16.3 



1.5 

51.5 

12.5 

49.0 

15.3 

39.5 

13.5 

27.3 

17.5 

25.0 

15.5 

21.8 

11.5 

Double 

1 

0.5 

18.0 

2.3 

11.3 

0.5 

11.0 

0.8 

11.5 

3.3 

10.5 

3.8 

10.3 

4.0 

-exp. 


1 

27.8 

2.8 

27.8 

2.0 

20.0 

2.0 

16.0 

4.0 

17.0 

4.8 

14.3 

4.0 



1.5 

34.0 

2.5 

31.8 

3.8 

30.8 

1.3 

17.0 

1.5 

18.3 

2.8 

17.5 

4.0 


5 

0.5 

29.0 

6.5 

20.0 

5.5 

19.0 

4.8 

17.8 

10.3 

17.3 

10.8 

16.8 

8.5 



1 

35.5 

7.3 

36.5 

7.3 

25.5 

4.3 

20.8 

12.0 

23.5 

13.5 

20.3 

14.8 



1.5 

42.3 

9.0 

38.8 

7.5 

37.3 

4.8 

23.8 

8.0 

25.0 

9.5 

23.0 

12.8 


10 

0.5 

36.0 

11.0 

26.0 

11.8 

22.5 

8.8 

23.0 

16.5 

22.3 

15.5 

19.8 

15.8 



1 

43.8 

13.5 

40.5 

13.3 

29.8 

10.3 

26.0 

17.8 

28.3 

23.0 

23.5 

20.0 



1.5 

46.5 

14.3 

42.5 

14.5 

40.5 

12.3 

29.0 

15.3 

29.5 

16.5 

27.5 

20.5 


and S = Id. The corresponding models are now denoted by ARi, * = 4,5,6 and due to 
the diagonal shape, we are dealing with two independent time series. In the second case, 
we generate realizations from the moving average model in (4.2) with B chosen from 


B4 


/O .5 0.5\ 

1 ^ 0.5 0.7)’ 


B5 


/O .5 0.5\ 

1 ^ 0.5 0.8 J’ 



0.5 

0.9 


and S = S 2 defined in (4.3) denoted by MAi, i = 4,5,6. Due to non-diagonal shape, 
we are dealing with two dependent time series in this case. To investigate the power 
performance of the tests, we consider the same non-linear time series models as above. 
First, we use independent GARCH(1,1) models, where the first process {A't_i,t € Zj 
follows GARCH(1,1) equation in (4.4) with uj = 0.01, a = 0.1, & = 0.2 and the second 
process {At_ 2 ,t G 1 } is generated by the same model, where w = 0.01, a = 0.1, but b is 
chosen from 


64 = 0.3, 


&5 = 0.4, 


65 =0.5 













Table 3. Actual size of and ip^ for nominal size a G {1%, 5%, 10%}, sample size n € {50,100,200}, bandwidth h = c - hcv for 
c€ {0.5,1,1.5} and non-linear models GARCHi^ TARi and RCAi^ respectively 


to 

to 


Model 

a . 

i : 

n : 

c 

1 






2 






3 






50 


100 


200 


50 


100 


200 


50 


100 


200 


p7l 

Tn 

pn 

Tn 

pn 

Tn 

pn 

Tn 

pn 


pn 


pn 


pn 


pn 

Tn 

GARCHi 

1 

0.5 

13.0 

3.5 

11.8 

3.8 

14.0 

3.3 

11.0 

3.5 

15.0 

2.3 

14.5 

2.5 

12.8 

2.5 

14.8 

4.5 

16.8 

6.5 



1 

15.0 

2.8 

14.5 

2.8 

14.3 

3.5 

12.5 

2.5 

21.8 

4.5 

17.3 

3.5 

16.3 

4.0 

16.0 

3.8 

17.8 

4.3 



1.5 

17.3 

3.0 

16.0 

3.5 

21.5 

3.0 

22.3 

3.8 

18.8 

4.0 

22.0 

3.8 

22.3 

3.8 

21.5 

4.0 

20.5 

3.5 


5 

0.5 

20.5 

11.0 

17.5 

9.8 

22.8 

10.8 

19.5 

12.0 

24.0 

11.8 

22.3 

9.8 

20.0 

11.5 

23.0 

13.0 

23.5 

14.0 



1 

21.3 

12.0 

20.3 

10.3 

21.5 

9.0 

18.8 

8.3 

28.0 

14.3 

23.3 

12.5 

22.8 

9.5 

22.5 

11.0 

25.8 

11.3 



1.5 

25.8 

7.3 

21.3 

8.8 

26.5 

11.3 

29.8 

11.0 

24.8 

12.0 

26.8 

13.0 

25.8 

12.0 

28.0 

14.0 

27.0 

9.8 


10 

0.5 

24.0 

16.8 

22.5 

15.3 

27.3 

19.0 

26.5 

18.0 

28.3 

17.8 

28.8 

19.5 

27.5 

20.8 

28.0 

19.8 

29.3 

21.5 



1 

28.3 

18.5 

27.0 

16.5 

28.8 

17.3 

25.8 

14.3 

32.5 

22.5 

27.5 

18.3 

29.3 

17.5 

29.0 

18.5 

31.3 

19.0 



1.5 

29.8 

13.5 

25.5 

13.8 

32.3 

17.5 

33.8 

17.8 

29.8 

16.8 

30.5 

19.5 

30.3 

18.8 

33.5 

20.5 

31.8 

17.5 

TARi 

1 

0.5 

7.0 

1.3 

8.3 

1.8 

6.0 

1.3 

8.5 

1.5 

8.3 

0.8 

8.3 

1.0 

10.5 

1.8 

10.0 

1.8 

11.0 

3.8 



1 

11.5 

1.5 

10.3 

1.0 

11.0 

1.5 

11.8 

1.0 

11.0 

1.8 

11.8 

1.3 

13.3 

1.3 

12.0 

2.0 

13.0 

2.3 



1.5 

15.0 

1.0 

12.0 

1.5 

12.0 

0.5 

16.0 

2.8 

13.5 

2.3 

12.3 

2.3 

11.8 

0.8 

13.8 

2.0 

13.0 

2.5 


5 

0.5 

13.0 

5.3 

14.0 

5.3 

11.5 

5.0 

15.0 

6.0 

12.3 

6.0 

14.0 

5.3 

15.8 

6.3 

15.8 

8.5 

19.3 

10.8 



1 

16.3 

5.3 

15.8 

6.3 

15.8 

6.3 

17.0 

4.8 

17.5 

7.0 

17.3 

7.3 

21.3 

7.0 

17.5 

6.8 

21.0 

9.5 



1.5 

20.8 

6.3 

17.5 

4.5 

18.0 

5.3 

23.5 

8.5 

21.5 

6.8 

17.5 

5.3 

17.8 

5.0 

19.5 

5.8 

22.5 

6.8 


10 

0.5 

16.5 

8.8 

19.3 

13.3 

15.8 

10.3 

19.5 

13.3 

19.0 

11.5 

18.0 

12.0 

22.3 

10.8 

21.3 

14.8 

23.8 

15.3 



1 

19.8 

11.5 

19.8 

12.8 

20.5 

11.5 

21.3 

10.3 

21.5 

12.5 

21.8 

14.0 

27.8 

14.5 

22.5 

12.8 

26.8 

16.3 



1.5 

23.8 

10.8 

20.5 

8.8 

22.8 

9.5 

29.0 

14.0 

26.3 

12.0 

24.0 

11.8 

22.0 

9.3 

24.0 

10.3 

27.5 

11.0 

RCAi 

1 

0.5 

6.5 

2.5 

7.0 

1.3 

8.3 

0.8 

9.0 

1.8 

11.3 

2.3 

11.8 

2.5 

10.3 

3.0 

12.3 

4.0 

11.8 

2.8 



1 

12.3 

1.0 

11.0 

2.0 

12.0 

2.5 

10.8 

1.8 

13.5 

2.5 

13.0 

0.8 

15.0 

3.8 

12.8 

0.8 

12.3 

2.0 



1.5 

12.5 

2.5 

13.3 

1.8 

14.3 

1.5 

15.5 

4.0 

12.0 

1.5 

15.8 

2.0 

15.8 

0.8 

16.5 

1.8 

17.3 

2.3 


5 

0.5 

13.0 

6.5 

13.8 

7.5 

14.8 

5.5 

14.3 

6.5 

18.8 

9.3 

18.8 

8.3 

17.8 

9.0 

18.5 

11.3 

19.8 

10.8 



1 

18.0 

7.0 

17.3 

5.0 

19.0 

5.8 

16.3 

7.3 

21.0 

8.5 

20.0 

7.3 

20.3 

8.5 

19.0 

6.8 

18.8 

8.3 



1.5 

18.8 

6.8 

17.8 

4.8 

21.5 

6.3 

23.8 

8.0 

18.3 

4.0 

22.0 

7.8 

20.8 

6.8 

25.5 

8.3 

24.5 

7.3 


10 

0.5 

20.5 

12.3 

21.0 

13.3 

17.5 

11.0 

20.8 

14.0 

24.0 

15.8 

23.3 

15.0 

22.8 

14.3 

25.5 

18.8 

24.5 

16.3 



1 

22.3 

11.8 

21.5 

10.3 

23.0 

12.8 

21.8 

12.3 

26.0 

15.0 

25.0 

15.0 

25.8 

13.8 

23.8 

14.3 

24.3 

14.8 



1.5 

23.5 

12.3 

20.8 

8.3 

27.0 

13.3 

28.3 

14.5 

21.0 

9.8 

25.5 

13.8 

26.0 

12.8 

31.8 

14.8 

29.5 

15.5 
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Table 4. Actual size of ipn for nominal size a £ 10%}, sample size n £ 

{50,100,200,500,1000,2000}, bandwidth hcv and models ARz, MAs, GARCHz, TARz and 
RCAz, respectively 


n 


Model 

a 

50 

100 

200 

500 

1000 

2000 

ARz 

1 

27.0 

20.0 

21.8 

19.3 

18.5 

18.3 


5 

33.5 

27.3 

27.3 

25.8 

23.5 

24.0 


10 

38.0 

33.5 

32.8 

30.0 

28.5 

27.3 

MAz 

1 

8.3 

6.5 

5.8 

6.3 

5.8 

4.8 


5 

13.5 

12.0 

11.8 

11.0 

10.0 

9.3 


10 

17.8 

15.8 

16.3 

17.0 

14.0 

13.3 

GARCHz 

1 

16.3 

16.0 

17.8 

16.8 

13.3 

13.5 


5 

22.8 

22.5 

25.8 

26.3 

19.5 

20.8 


10 

29.3 

29.0 

31.3 

30.0 

27.3 

26.0 

TARz 

1 

13.3 

12.0 

13.0 

10.0 

8.5 

7.8 


5 

21.3 

17.5 

21.0 

15.5 

15.3 

14.8 


10 

27.8 

22.5 

26.8 

20.3 

21.3 

17.8 

RCAz 

1 

15.0 

12.8 

12.3 

13.5 

11.8 

11.5 


5 

20.3 

19.0 

18.8 

19.8 

18.8 

17.8 


10 

25.8 

23.8 

24.3 

25.3 

24.3 

22.0 


referred to as model GARCHi, i = 4,5,6. Two (centered) independent threshold AR 
models under the alternative are generated by equation (4.5), where {Xtp,t £ Z} cor¬ 
responds to coefficients (a(l),a(2)) = (—0.2,0.!) and those for {Xt^ 2 ,t £ Z} are chosen 
from 

a4 = (-0.3,0.2)'^, 05 = (-0.4,0.3)^, Og = (-0.5,0.4)"^. 

Remark in this context, that for the asymptotic test we would observe a larger 
absolute power since this test is very liberal in the considered situations. In the supple¬ 
mentary material (cf. Jentsch and Pauly [17]), we therefore present its power performance 
in comparison to all three randomization tests by a comparison of actually achieved size 
obtained from the simulations in Section 4.2 with rejection rates in Tables 5-6 instead 
of using the nominal size. There it can be seen that the power behavior of all three tests 
(measured as power in comparison to actual size) appear to be quite similar for all models 
under investigation. 

4.4. Discussion 

From Tables 1-4, it can be seen that the asymptotic test ipn has difficulties in keeping the 
prescribed level and tends to over rejects the null systematically for all small (n = 50) 
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Table 5. Power of for nominal size a € {1%, 5%, 10%}, sample size n G {50,100,200}, band¬ 
width h = c - /lev for c G {0.5,1,1.5} and autoregressive models AR4, AR5 and ARe with Gaus¬ 
sian, logistic and double-exponential distribution of the innovations 


{et} 

a 

n: 

c 

ARi 



AR 5 



AR% 



50 

100 

200 

50 

100 

200 

50 

100 

200 

Gauss 

1 

0.5 

1.5 

3.3 

10.8 

5.3 

9.0 

35.0 

4.8 

19.3 

64.3 



1 

3.0 

5.3 

10.0 

4.5 

11.8 

37.0 

6.8 

19.0 

67.0 



1.5 

3.8 

2.5 

11.5 

3.8 

10.0 

36.8 

5.3 

20.5 

67.5 


5 

0.5 

7.5 

14.3 

32.3 

16.5 

29.0 

64.8 

19.0 

49.5 

89.0 



1 

8.5 

15.3 

25.8 

14.5 

33.8 

70.3 

23.3 

51.5 

91.0 



1.5 

12.5 

12.3 

33.0 

12.8 

29.5 

71.0 

20.3 

48.3 

89.3 


10 

0.5 

16.5 

24.0 

46.3 

28.5 

46.8 

80.8 

31.5 

72.5 

95.0 



1 

17.0 

26.0 

39.0 

26.0 

47.5 

83.8 

39.5 

71.8 

96.3 



1.5 

18.8 

25.0 

46.5 

23.3 

43.3 

84.5 

36.0 

64.3 

97.3 

Logistic 

1 

0.5 

3.3 

5.0 

10.8 

6.3 

10.0 

36.8 

8.0 

24.0 

64.5 



1 

3.8 

4.8 

11.0 

6.5 

12.5 

31.5 

7.3 

16.3 

70.0 



1.5 

2.8 

5.8 

10.5 

6.0 

12.3 

34.3 

10.3 

22.3 

66.5 


5 

0.5 

10.0 

17.0 

29.3 

16.5 

31.5 

65.5 

22.8 

55.5 

88.3 



1 

11.3 

16.5 

26.5 

20.3 

33.3 

64.5 

21.8 

53.0 

90.0 



1.5 

9.0 

18.3 

26.8 

16.3 

30.0 

64.5 

24.5 

52.5 

91.0 


10 

0.5 

18.0 

26.5 

45.3 

25.8 

48.5 

77.3 

37.8 

69.8 

94.8 



1 

18.5 

26.0 

44.0 

29.5 

50.5 

81.0 

35.5 

71.5 

96.0 



1.5 

18.3 

26.5 

39.0 

26.5 

44.0 

81.5 

39.5 

70.5 

98.3 

Double 

1 

0.5 

4.0 

5.0 

9.8 

6.3 

10.3 

39.8 

10.3 

24.0 

59.3 

-exp. 


1 

6.0 

6.0 

8.5 

5.8 

13.8 

38.3 

00 

bo 

23.5 

65.5 



1.5 

4.5 

6.8 

12.5 

8.0 

15.0 

35.3 

9.3 

23.5 

69.0 


5 

0.5 

13.3 

16.5 

27.0 

18.5 

31.5 

67.3 

27.3 

54.8 

86.8 



1 

14.5 

16.8 

24.0 

17.3 

34.3 

66.8 

23.8 

55.3 

91.0 



1.5 

12.0 

17.0 

32.8 

19.0 

38.3 

66.5 

24.0 

50.3 

89.0 


10 

0.5 

22.3 

25.5 

43.0 

28.8 

48.3 

82.0 

41.0 

75.0 

93.8 



1 

23.3 

28.5 

40.5 

31.0 

49.3 

82.5 

39.0 

71.8 

95.8 



1.5 

21.0 

23.3 

43.8 

31.0 

51.8 

79.5 

39.0 

66.5 

97.0 


and moderate (n = 100) sample sizes. Its performance is not even desirable for larger 
sample sizes (n > 200). This is the case for all linear and non-linear as well as all three 
innovation distributions under consideration. In Table 1, especially for the most critical 
autoregressive models AR 3 , where the AR coefficient is near to unity and the correspond¬ 
ing spectral densities have a non-flat shape, the null approximation of ipn is extremely 
poor and the performance is unacceptable (see Table 1, right panel). Nevertheless, this 
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Table 6. Power of for nominal size a £ {1%, 5%, 10%}, sample size n £ {50,100,200}, band¬ 
width h = c - hcv for c £ (0.5,1,1.5} and moving-average models M% 4 , MA^ and MAq with 
Gaussian, logistic and double-exponential distribution of the innovations 


{et} 

a 

n: 

c 

MAi 



Mds 



MAe 



50 

100 

200 

50 

100 

200 

50 

100 

200 

Gauss 

1 

0.5 

3.3 

7.0 

14.3 

4.5 

12.3 

34.8 

12.8 

29.8 

65.5 



1 

4.8 

11.0 

17.8 

7.8 

19.5 

38.0 

12.3 

38.3 

69.0 



1.5 

3.3 

9.5 

20.8 

9.5 

19.0 

44.3 

14.0 

34.8 

76.5 


5 

0.5 

10.0 

17.0 

31.0 

13.8 

29.3 

54.0 

31.5 

50.0 

84.8 



1 

12.8 

21.5 

34.3 

21.8 

40.3 

62.8 

29.8 

62.0 

87.0 



1.5 

12.5 

20.8 

39.0 

23.3 

39.0 

66.3 

29.3 

62.0 

90.0 


10 

0.5 

19.8 

26.8 

41.3 

22.8 

41.0 

64.5 

41.5 

63.8 

91.5 



1 

20.0 

34.0 

47.3 

31.8 

53.0 

74.5 

43.0 

72.0 

94.3 



1.5 

23.3 

30.5 

50.0 

35.3 

52.8 

78.8 

43.3 

74.5 

95.5 

Logistic 

1 

0.5 

2.5 

3.8 

8.3 

3.0 

10.8 

16.0 

9.3 

14.8 

30.0 



1 

2.5 

7.0 

12.3 

5.3 

11.8 

14.5 

10.8 

20.0 

36.8 



1.5 

2.5 

5.0 

9.0 

5.0 

11.0 

21.0 

10.3 

14.8 

40.0 


5 

0.5 

13.5 

13.3 

19.0 

13.5 

23.5 

29.5 

22.5 

29.0 

49.0 



1 

11.8 

13.3 

22.8 

14.8 

26.0 

34.3 

24.3 

39.3 

59.5 



1.5 

10.8 

12.3 

22.0 

15.3 

24.8 

39.3 

21.0 

33.5 

61.8 


10 

0.5 

19.0 

21.3 

28.3 

21.3 

34.3 

41.8 

31.3 

39.3 

59.8 



1 

19.8 

22.5 

33.0 

24.0 

36.8 

47.3 

32.3 

51.0 

71.3 



1.5 

18.5 

19.8 

31.0 

21.5 

37.0 

51.5 

34.5 

45.5 

72.0 

Double 

1 

0.5 

3.5 

4.8 

9.0 

6.3 

11.8 

18.5 

11.3 

19.8 

35.0 

-exp. 


1 

6.5 

10.5 

11.3 

8.0 

13.3 

23.0 

10.3 

20.5 

37.3 



1.5 

5.3 

9.8 

13.5 

7.0 

12.0 

22.3 

12.3 

23.5 

43.5 


5 

0.5 

12.5 

14.8 

20.5 

19.5 

24.5 

32.8 

26.8 

34.8 

51.8 



1 

13.3 

19.8 

25.3 

19.5 

26.0 

39.3 

23.8 

39.3 

55.5 



1.5 

18.0 

21.5 

24.8 

18.3 

25.8 

40.8 

25.0 

40.5 

64.5 


10 

0.5 

21.5 

23.8 

28.3 

29.8 

36.3 

42.3 

35.8 

45.0 

60.5 



1 

23.0 

27.3 

36.3 

29.3 

35.8 

48.5 

33.0 

50.3 

64.5 



1.5 

26.5 

34.0 

36.5 

26.8 

36.5 

49.3 

34.0 

50.3 

73.8 


generally poor performance is not surprising since the slow convergence speed of L 2 -type 
statistics is already known, see for instance Paparoditis [21], and Table 4 in this paper. 

In comparison to that, the randomization test performs better than ipn- For all 
models under the null and innovation distributions under investigation, the randomiza¬ 
tion test holds the prescribed level more satisfactorily than ipn- 

Especially for Gaussian innovations and linear time series, the usage of can be 
recommended for all dependent and independent settings and sample sizes. In these 
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cases, a close inspection of Tables 1-3 shows also that the bandwidth choice seems to 
have only a slight effect on the behavior of the randomization test, where this choice 
appears to be more crucial for To demonstrate this, compare for instance in Tables 1 
and 2 the performance of model AR^ for n = 200 and a = 10%, where the range of 
the actual size is from 22.5% to 40.5% for tpn and from 8.8% to 12.3% for over the 
bandwidths 0.5hcv,hcv and 1.5/icv- 

For linear time series with non-Gaussian innovations the performance of ipn is still 
worse than that of (/?*. However, here the randomization test is slightly more effected by 
the choice of the bandwidth and its accuracy under Hq often needs larger sample sizes 
(n > 100 or even n > 200). 

For non-linear time series, see Table 3, a similar observation can be made. The perfor¬ 
mance of ipn is poor for all sample sizes from 50 to 200 and again keeps the prescribed 
level much better. In particular, for the TAR models, the control of the nominal size is 
quite accurate. For the RCA and GARCH models under consideration, which are not 
covered by our Assumption 2.1, the performance of is still quite good and, in partic¬ 
ular, improved in comparison to (fn- However, the finite sample performance seems to be 
more effected by the bandwidth choice for RCA models and still tends to overreject 
the null for all GARCH models. 

Tables 5-7 show the power behavior of the randomization test , where we compare 
its power to its nominal size. When studying the panels with increasing sample sizes 
(from left to right) the consistency results of Theorem 2.2 and Corollary 3.1 under the 
alternative can be confirmed by the simulations. In particular, for the non-linear RCA 
and GARCH models a typical consistency behavior can also be observed. Similar to the 
situation under the null, for Gaussian innovations and linear time series the bandwidth 
choice does not seem to be crucial for but for other innovations and non-linear 
time series it has a considerable effect on its power behavior. To this end, we think 
that the typically applied cross-validation selector leads to quite adequate finite sample 
performances. 

Our simulation experience may be summarized as follows: 

• The randomization technique makes sure that ip’^ keeps the prescribed level well for 
very small sample sizes (especially for most linear time series and TAR models under 
consideration) as shown in Tables 1-3. Its performance is in general significantly 
better than those of the asymptotic benchmark test (p„. 

• Even for larger sample sizes up to n = 2000 the asymptotic test (pn cannot be rec¬ 
ommended due to its very slow convergence as emphasized in Table 4. 

• In comparison to other resampling methods applied in time series analysis, the ran¬ 
domization technique used here has the big advantage that their performance does 
not depend on the choice of any tuning parameter in addition to the bandwidth. 
The choice of the bandwidth can be done automatically by standard methods and 
does not seem to be as crucial as for pn- 

• The performance of (^* becomes even more excellent if we compare its behavior for 
the very small sample size of n = 50 with the poor performance of the unconditional 
test. 


Testing equality of spectral densities 


27 


Table 7. Power of for nominal size a £ {1%, 5%, 10%}, sample size n £ {50,100,200}, band¬ 
width h = c ■ hcv for c £ {0.5,1,1.5} and non-linear models GARCHi, TARi and RCAi for 
i = 4,5,6, respectively 


Model 

a 

i: 

n: 

c 

4 



5 



6 



50 

100 

200 

50 

100 

200 

50 

100 

200 

GARCHi 

1 

0.5 

5.3 

6.5 

8.5 

12.3 

22.8 

42.0 

21.0 

53.5 

86.3 



1 

4.3 

7.5 

11.8 

12.3 

23.8 

46.5 

26.5 

56.0 

83.0 



1.5 

4.0 

5.3 

11.3 

12.8 

24.0 

46.8 

25.5 

55.3 

83.5 


5 

0.5 

12.8 

16.3 

20.5 

26.3 

35.8 

61.3 

41.5 

71.3 

95.0 



1 

10.3 

17.8 

25.3 

24.5 

38.5 

66.8 

47.3 

72.3 

92.8 



1.5 

12.8 

13.8 

21.0 

28.3 

40.5 

66.3 

49.5 

73.5 

94.3 


10 

0.5 

19.8 

24.5 

27.3 

37.0 

48.8 

70.3 

56.5 

82.5 

96.3 



1 

16.8 

25.3 

34.8 

32.8 

50.0 

77.3 

58.5 

81.0 

95.8 



1.5 

21.3 

23.0 

33.0 

37.3 

51.5 

74.8 

59.0 

81.3 

96.8 

TARi 

1 

0.5 

2.5 

2.3 

3.5 

2.0 

1.3 

2.8 

2.0 

3.8 

6.8 



1 

1.5 

1.3 

1.8 

1.3 

0.8 

3.8 

2.0 

3.5 

5.8 



1.5 

1.5 

1.5 

2.8 

2.0 

1.3 

1.0 

1.8 

1.5 

3.3 


5 

0.5 

6.3 

8.3 

8.0 

8.0 

6.0 

12.8 

6.5 

11.3 

21.3 



1 

6.5 

5.5 

4.5 

6.3 

4.8 

9.8 

6.3 

10.5 

18.5 



1.5 

5.8 

5.5 

7.8 

6.3 

5.5 

5.8 

5.8 

6.5 

12.5 


10 

0.5 

13.3 

12.5 

13.3 

12.8 

10.3 

22.3 

14.0 

22.5 

30.8 



1 

12.8 

10.8 

10.5 

11.5 

9.8 

15.8 

11.8 

17.5 

27.5 



1.5 

11.5 

9.3 

12.5 

11.8 

11.0 

9.0 

13.8 

14.0 

19.3 

RCAi 

1 

0.5 

1.8 

1.3 

2.3 

2.5 

3.0 

2.8 

2.3 

6.0 

10.5 



1 

1.3 

2.5 

2.3 

1.5 

3.8 

2.5 

2.8 

5.8 

12.0 



1.5 

2.3 

0.8 

1.3 

1.5 

1.5 

2.5 

3.8 

5.8 

12.0 


5 

0.5 

6.0 

7.3 

7.0 

6.8 

11.3 

10.8 

8.5 

15.0 

20.3 



1 

5.8 

7.0 

8.3 

6.5 

11.3 

9.3 

13.0 

15.3 

27.3 



1.5 

6.8 

7.5 

5.5 

6.3 

8.8 

9.8 

10.3 

16.0 

24.5 


10 

0.5 

14.8 

12.5 

12.5 

12.3 

14.8 

16.5 

17.3 

23.3 

30.3 



1 

10.0 

14.3 

12.3 

12.5 

16.5 

15.0 

23.3 

23.8 

39.3 



1.5 

14.0 

13.0 

13.5 

11.0 

14.8 

17.0 

19.8 

22.8 

33.0 


• In comparison to that the only other known and mathematically analyzed resampling 
test for iJo, the bootstrap test of Dette and Paparoditis [9], needs sample sizes of 
n > 512 for gaining comparable results (see Table 1 in their paper). Remark that 
they have also modeled a VAR(l) model under the null with Gaussian innovations. 
When comparing their results with ours (see Table 1 in this paper) note that p = 0 
in their paper corresponds to our model AR^, but with the (slightly less critical) 
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choice 



• Furthermore, the power performance of improves with increasing sample size (as 
usual under consistency). The simulations in the supplementary material even show 
that the power behavior (measured as power in comparison to actual size) is similar 
to that of the benchmark test ipn and all other randomization tests. 

• Finally, to sum up, the randomization procedure helps to hold the prescribed level 
under the null more satisfactorily and does not forfeit power under the alternative 
in comparison to the unconditional case. 

5. Final remarks and outlook 

In this paper, we have introduced novel randomization-type tests for comparing spec¬ 
tral density matrices. Their theoretical properties have been analyzed in detail and we 
have also studied their finite sample performance in extensive simulation studies. The 
asymptotic behavior under the null as well as for fixed and local alternatives have been 
developed for non-linear time series under a joint cumulant condition (2.1), whereas for 
linear processes finite 32nd moments are required. Although these conditions seem to 
be rather strong, we like to mention that only a few mathematical results on this topic 
can be found in the current literature, where most of them are only developed for the 
one-dimensional case. Note, that such, or stronger linearity or even Gaussianity condi¬ 
tions are typical assumptions in recent time series publications, see, for example, Eichler 
[11], Dette and Paparoditis [9], Dette and Hildebrandt [8], Jentsch [15], Jentsch and 
Pauly [16], Preuss and Hildebrandt [25] or Preuss et al. [26]. Moreover, one of the main 
contributions of the current paper is the theoretical justification of the randomization 
approach, which has yet not been achieved in the time series context. 

In future research, we aim to get rid of the last tuning parameter (i.e., the bandwidth) 
by studying a different test statistic that is based on integrated periodograms rather than 
kernel spectral density estimators. In doing so, we also plan to relax the conditions on 
the process in the non-linear case and to investigate to what extend our theory can still 
be established. As pointed out by one referee, a promising approach may be given by 
substituting our joint cumulant condition by physical-dependence-type-conditions intro¬ 
duced in Wu [32] and further studied in for example, Liu and Wu [19] or Xiao and Wu 
[33]. However, to our knowledge the theory is currently not available in full strength for 
our setting, that is, for multivariate fc-sample problems with dependencies and triangular 
arrays. The latter is required to derive the local power behavior of our test statistic. 

6. Proofs 

In the following and for a better Incidity of the proofs, we will use the abbreviate notation 

^w,mi ,m2'-,k — ^pw—mi ,pw — 7n2 ) ^od fwi,W2,'mi,m2;k — fpwi—mi ,p'W2 — m2 i.^k ) • 
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Proof of Theorem 2.1. Under the Assumptions 2.1 and 2.2, which accord with As¬ 
sumptions 3.1 and 3.3 (with taper function h(t) = 1) in Eichler [11], Theorem 3.5 in his 
paper is applicable to the test statistic r„, where the also required Assumption 3.2 is 
obviously satisfied due to the quite simple shape of r„. Also, by direct computation, we 
get 


ul/2 nn P-1 

E{Tr,) = — ^ 

mi ,7712=0 ki ,/C2 


^^2 


J2 = l \l-=l 


X Eil,, 

,mi,m2]kilj2,mi,m2;k2) 


nh?!^ 


P-1 

E E^^ 

7712=0 k 


2 U-UJk 




ll j'2 = l 

which is asymptotically equivalent to 




For the variance, we have 
Var(T„) 


n 


%r r 


f—nJ—n 


p-1 

E E Kh{uJ - uJkf)Kh{aj - LOk^)Kh{X - uJk3)Kh{X - ujkj 

mi ,7712,7713,7714 = 0 fci ,/C2,fc3,fe4 


1 


y f y (1 

'\ll >12.13 04 = 1 \l'l=l / 


X j y ] (1 Q<^l3r2)(l 


^T’2 = l 


X {E{lji,mi,m2-,kilj2,m,i,m2;k2lj3,m3,7niik3lj^^m3,m4-,k4) 


-E{In ,7711,7712 ;fci ^j2 ,7711 ,7712 ;fc2 ) 


^ ^ ,7713,7714 ;fc3 '^^4,7713 ,7714 (^4 )} dA. 
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In evaluation of the difference of the expectations above, the only asymptotically non¬ 
vanishing cases are uJki = Wfcj Wfcs = = Wfc 4 ^ ^k-^ = Wfcj, = -ujk^ ^ LUk^ = 

—ujki and ujk-i = —i^k^ Wfc 2 = ~^fc 3 , where all of them make the same contribution, which 
enables us to consider the first case only equipped with a factor 4. In this situation, we 
have 




4;fci /: 




■,kif: 


32,34:,mi,rrii 


fc 2 fr. 


32,34,m 2 ,m 3 


;k2 ^kiks 5k2k4 j 


which results in 


Var(r„) 

^ /' \ ^ I 

= — ^ / Khiuj-uJki)Khiuj-ujk3)duj\ — 

^ ki,k2^-^-^ ^ ^ 


jl ,12 03 04 = 1 


I f 7-^ 

X 1 / ^ 3 31,33,mi,m4-,kiJj2,j4,mi,m4-,k2 I 

\mi,m4=0 / 

/ p-i _ \ 

^ ( 'y ^ I3l,33,m2,m3-,kifj2,34,m2,m3-,k2 1 +o(l)- 

\m2,m3=0 / 


The term on the right-hand side is asymptotically equivalent to 


Tn =Bk 


j (-l + '?^Jii2)(-l + 9<5j3j4)|tr(Fjjj3(u;)F^,j^(a;)^)n dw, 


il ,32,33,34 = 1 

which concludes this proof. 


□ 


Proof of Theorem 2.2. Note Hrst that the -x/uTi-consistency of the kernel density 
estimators, see, for example. Equation (5) in Eichler [11], implies p, — iiq = op{h^^‘^). 
Thus Slutsky’s lemma together with Theorem 2.1 show that under Hq 


- h-V2 


^3 V 


T-O 


iV(0,l) 


which corresponds to the limit result for Qt on page 976 in Eichler [11]. Together with 
the consistency of this proves the asymptotic exactness of (/?„ under the null hypoth¬ 
esis. In contrast to that, similar to Eichler’s Theorem 5.1, it can be verihed by direct 
computations and Assumption 2.2 that converges in probability to -l-oo under the 

alternative Hi. Thus, the same hold true for Tn — which implicates consistency 

of the test. □ 
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Proof of Theorem 2.3. Remark that the imposed uniform mixing conditions still 
imply the i/uTi-consistency of the kernel density estimators for the processes. Hence, the 
estimators for mean and variance can be replaced by their limit values and we can apply 
Theorem 5.4 in Eichler [11] with the special function 


vk(f)= (vec(F„-i^F,y) ) 

\ \ ^1=1 / / l<r<g 

Noting that 

= vec((l{(r,s) = - -vec((l{(r,s) e 


( {{r,s) eE: 3c with {r,s) = {cp + i,cp +j)}, for(j,j)eE;, 

1 ) ~ 1 

0, elsewise 

and with S = ljg_]^{(r, s): 1 + (c — l)p <r,s< cp} that theorem implies, that our stan¬ 
dardized test statistic converges in distribution to a normally distributed random variable 
with mean v and variance Tq . Hence, the results follows. □ 


a4'(Z) 

dZ,. 


6.2. Proofs of Section 3 


Throughout this section, E*(-) := ... ,X„), Var*(-) := Var(-|X^,... ,X„) and 

Cov*(-,-) := Covf-, - lx,...., X„) denote mean, variance and covariance conditioned on 
the data, respectively. 


Lemma 6.1 (Brillinger [3], Theorem 4.3.2). Suppose that the mixing condition (2.1) 
is satisfied for all k < s Then for any k £ {1,..., s} and wi,..., G [0, 27 t], we have 


CUm( Jai (Wi ),..., Jafc (Wfc ) ) 


„fc/2 Jk 


(Wl, 




j = i 



( 6 . 1 ) 


where Jai<-v) denotes the ath component of the d-variate discrete Fourier transform Jfoj), 
see (1-4), is defined as 


(27t) 


k-1 


E 


,{ti,..., tk-i) exp(i(tiWi H-h tk-iUJk-i)) 


— l — OO 


the kth order cumulant spectral density and fk is the (dx dx ■ ■ ■ x d) (n times) kth order 
cumulant spectra of the d-variate time series (X^.t £ Z). 
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Lemma 6.2. It holds 

9 9 

(^) ^ ~ .'"■fci (r))(l ~ 1^32,Trk 2 (r)) — ~Q P Q ^ ^ {r)^j 2 ,'^k 2 

r—1 r—1 

Q 

(il) ((^ ~ ^^ii,7rfcj(r))(l ~ Q''5j2,'7rfc2C’^))) “ “I" ^ ^iii2 )'^| fci |, 1^2 h 

r—1 

Q 

(iii) ^ (l’l))(l 9'^j2,7I-fc2(ri))(l 9'^j3.'n’fc3(r2))(l .7rfc4 (l’2))) 

ri ,r2 = l 

= q (—1 + 9'^ili2)(~l + Q^j3j4)^\ki\,\k2\^\k3\,\k4,\ (1 ~ ^IfclMfesl) 

+ 9 ^ _ 2^ (1 ~ '^iii3)(l ~ ^hoS^ 

X <^|felMfe3l<^|fc2|.|fc4|(l-'^|fel|.|fc2|) 

+ 9 ^^hii^hh + ^ (1 ~ '^iii4)(l ~ ^1233^ 

X ^|felMfe4|'^|fc2|.|fc3|(l“^|fel|.|fc2|) 

+ C'(5|fci|_|fe2Mfe3Mfc4h 


/or some constant C < oo, where ^|fci|,1^21.1^31,1^41 = 1 */ l^il = l^al = l^al = 1^41 and 0 
otherwise. 


Proof. The first assertion (i) follows from Y11-=i^hTTk(r) = 1- Due to P(7rfc(r) = j) = 
for all k,r and j, we get 


^ i^3l,'^kir)Sj2,Trk(r)) 

^ (^3l,'^kiiri)^j2,Trk2ir2)) ~ \ ^ (*^91 .’^fc (’’l) ^t2 .TTfc (’’2) ) i 


\ki\ = \k2\,ri =r2, 

|fci| = \k2\,ri ^r2, (6.2) 
i^3i,^kAr))^* (^32,^k2ir)), I fell 7^ 1 ^ 2 ! 


-5, 


1 


= 


l^il = 1^2!,?'! =r 2 , 

i'^-^ 3132 )^ l^il = |^2|,?'i 7^r2, 

l^il 7 ^ 1 ^ 2 ! 


qiq-iy^ 

1 

7^ ’ 


which yields (ii). Now consider (iii). First, set Zj^k,r := 1 — 9^j,7rfc(r) and note that 
E*{Zj^k,r) = 0 for all j, k and r. Therefore, we have to consider the following non-vanishing 
cases (a) |fci| = \k 2 \ ^ Ifca] = \k4\, (b) |A:i| = l^al ^ jfel = \k4\, (c) \ki\ = \k4\ ^ |/c2| = l^al 


Ih-i 
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and (d) |fci| = |A: 2 | = l^al = |fc 4 |. The first case (a) gives 

q^(-l + qSj^j2){-l + 9<^i3j4)<^|fcl|,|fe2|'^|fe3l.|fc4|(l “ '^Ifcll.lfesl) 
which follows immediately from (ii). Case (b) becomes 

Q Q 


E 


^ i^h,ki,riZj2,k2,riZj^^k3,r2^j4,ki,r2) — ^ ^ ^ (-^Ji .fci ,ri,fei ,r2 .^2 .»'l ^^4 .fc2,i’2 ) 


i’l,r -2 = l r'l,r '2 = l 

and together with ( 6 . 2 ), we get 


^ iZji,ki,riZj2^ki,r2) — 1 + 9<^rir2<^ji 43 "I" ■ 




and an analogue result holds for ^I*(^j 2 ,fe 2 , 1 - 1 ^ 44 ,fc 2 .r- 2 )- 
up, the contribution of case (b) becomes 


q ( -1 + g(^4i43<5j2j4 + 


9-1 


(1 


4143 


By multiplication and summing 
)(1 — '^4244)^ 


and case (c) contributes analogously. □ 

With these results, we can analyze the conditional expectation and variance of T*. 

Theorem 6.1 (Conditional mean and variance of T*). Under the assumptions of 
Theorem 3.1, it holds 

i?*( 7 :)=h-i/v:+op(i), 


where 


^ l-Tt P-1 ^ 9 

Tn'-~~ I 'y ] 'y ] ~ ^k)— 'y ] (—l+9'^4i42)-^P4l-"U.P4l-"*2(l^fc) 


mi,7712=0 k 


iio'2=i 


^ ,p_72 — 


(6.3) 


Moreover, we have 

Var*(T:) = ^ E (rKniu 

kuk2 

\ki\^\k2\ 


UJki)Kh{u} 


Ujk2 ) dw 


2 


P-1 1 / 

E ^ E ( — 1 + 9 <^ 4143'^4244 

mi,m2,m3,m4=0 ^ 4l J2.43.i4 = l ' 


X 
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X L 


ji ,mi ,m2',ki Ij2 ,7711,m2',k2 


^ ^j3,m2,7n4,',k\ij4,,7n2,mi',k2 


=:Tf+Op(l). 

Proof. By using Lemma 6.2(ii), we get 

/ll/2 


^71 P 1 1 

E*{T:) = — ^ 

^ 0 fci,fc2 ^ 


^ Z! ^*( .’^'=1 (’■))(! -9^i2.7r,,(r-)) 

ji,i 2 =i \r—i y 


^ ,?71i ,77l2;fci-^_72 ,7712 ;A:2 


/ll/2 


P-1 


y^Kh{u} - u}ki)Khiu} - LUk^)- 

' mi,7712=0 ki,k2 


^ Z! ( 1 + 9^iii2)'^|fci|,|fe2| 
tl J2 = l 


^ ^ji,m7,m2',ki^j2,mi,77i2;k2 


Furthermore, the case fci = —^2 is asymptotically negligible, which yields the first asser¬ 
tion of Theorem 6.1. Considering the conditional variance of T*, we get 


Var*(rf) 

_ h r 

J-n 


P-1 


E E Kh{uj - uJk^)Kh{uJ 

'■ mi ,m2,m3,m4=0 k\,k2,ks,k4 


^k2 )^h 


ujk3)Kh{X-uJk4) 



X 


E 


q 

^ ' COV •^j2,fe2,I’l I 


ll J2 J3 J4 = l ri,r2 = l 


,k3,7'2^3i,ki,r2) 


^ ^ji7'fT^\-,'^2'-,ki^j2-,Tn\,m2\k2 


^ ^j3,'f^3,fT^4',k:i^j4,m3,m4^-,k4 dw dA. 
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By using Lemma 6.2(ii) and (iii), only the cases |fci| = l/ca] ^ \k 2 \ = \k 4 \ and |fci| = |/c 4 | 
1 ^ 2 ! = l^al play a role asymptotically. More precisely, only the four cases fci = fcs ^ ^2 = ^4 
and fci = —fca ^ ^2 = —ki and fci = ^4 ^ ^2 = fca and ki = —ki ^ ^2 = —^3 do not vanish, 
but all of them make the same contribution asymptotically. This gives a factor 4 and we 
get up to an op(l)-term 

Var*(rf) 

Tt P-1 

E E Kh{uJ - uJki)Kh{uj - LOk^)Kh{X - u}ki)Kh{X - 

mi,m 2 ,m 3 , 1714—0 ki,k 2 

\ki\^\k2\ 

I 9 / 

E/ ( 

il J2 J3J4 = 1 ^ 

+ ^ (1 ~ '^4143 )(1 ~ 

^ ^jl,mi,m2;ki^j2,mi,m2',k2 

^ ^33,m3,m4',ki^j4,m3,m4\k2 dw dA, 
which concludes this proof. □ 

Theorem 6.2 (On E{yd^), Var(/r*), and Var(r*^)). Under the assumptions of 

Theorem 3.1, it holds 



) =/U* + o(h^/^) and Var(^*) = o(/i), 
E{t*^) = T*^ + o{l) and Var(T*^) = o(l) 


(6.4) 


as n —?► 00 . 


Proof. Because arguments are completely analogue, we prove only the more complicated 
part that deals with r*^. First, by introducing the notation Jj,m-,k = St=i , 

we get 


Ah _ / \' 

E(,t*'^) = — / Kh{u}-u}ki)Kh{uj-u}k2)duj] 

” h, \J-rc / 


ki ,/C2 


p-1 

E 


mi,m2,m3,m4=0 41,42,43,44 = 1 


q" 


E (-1 + '?'^ 4 i 43 ^: 


J2J4: 


q-1' 


(1 '^4143 )(1 '^4244) 
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(6.5) 



The last expectation above can be expressed in terms of cumulants and we have 


EiJi • • • Jg) = En cum(Jii i G B), 


( 6 . 6 ) 


cr B^<7 


where we identify the DFTs with Ji,..., Jg for notational convenience, cum( i G B) \s 
the joint cumulant of [Jg iG B), a runs through the list of all partitions of {1,..., 8} and 
B runs through the list of all blocks of the partition a. Note that the largest contribution 
in (6.6) is made by products of four cumulants of second order. Moreover, all of those 
combinations where cumulants cum( J...;fc., J.^.jfe^) or cum(J.^.;fe^, J._.;fc^.) for i^j occur, 
contain a factor Sk^k 2 oi 5ki-k2i respectively, but these cases are excluded in the sum¬ 
mation in (6.5). The cases with cum( or cum( contain a factor 

which causes the sum over ki to collapse and, therefore, they are of lower order and 
asymptotically negligible. The only two remaining products consisting of four cumulants 
of second order are by (6.1) 


cum( 




J3 




and 




Now, by taking the sums over mi,m 2 ,mg ,7714 of both expressions, we get 



(6.7) 


and 
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respectively. Asymptotically equivalent to (6.5), this results in 


Bk 




X {tr(Fjjjj (wfcj )Fj2(wfcj) ) tr(Fj3j3 (wfej (wfej) ) 


+ |tr(Fj3j3(a;)Fj2j^(w))! }da;, 

which shows the first assertion of (6.4). For its second part, we have 
Var(r*2) = ^^ ^ ( f Kh{uj - uJkjKhiuj - duj') 


^1,^2 k3,k4 

\ki\^\k2\ 1^31/1^41 


J Kh{X-uJk 3 )Kh{X-uiki)duj^ 

p-i 1 / 

y~! 71 ^ [~^ + qSj^j^Sj. 


mi,...,ms=0 ji,...Js = i 


+ '^jlj3)(l ^j2j4) 

^ + d^jsjY^jeh 


+ ^_2^(1 ^t5j7)(l Sjgjg) 




^ Ji'i J' 


J3,m4;fei ’Jj2,m,i;k2Jj2,m2\k2 


X '^i4.m3;fe2'^i4,"l4;fe2'^j5,m5;fc3 
^ '^i5."J6;fe3'^i7,"l7;fe3 (6-8) 


^ Ji-T .m.s'.k'i J id .m.^'.k.A J ’ 


J7 .mg;fe3 >7.ms;^4 .7«6;fc4 


^ Jjs,m7;k4Jjs,ms-,k4) 


B{Jji,mi-,ki Jji,m2-M '^i3,”i3;fel 


^ <^^3 ,m4;fci ,mi ;fc2 
^ '^J2 ,7J72;fc2'^j4,"i3;fc2 
^ Jj4,m4-k2) 


X E{Jj^^m5-,k3Jj5,me-,k3 
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X Jj7,mr-,k3JjY,ms;k3 
X Jje,m5-,k4Jje,Tne\ki 
X Jj3,mT,kiJig,mg-,ki)} 

and similar to the computations for we are able to express the 16th and both 

8th moments with cumulants and the difference above becomes 

{E(Ji ■ ■ ■ Jie) - E(Ji ■ ■ ■ Ji)E{J^ ■ ■ ■ Jie)} = En cum(Ji: iGB), (6-9) 

where we identify the DFTs Ji^rn\k from above with Ji,..., Jiq for notational convenience 
and now a ruirs through the list of all partitions of 16} that can not be written 

as cumulants depending on subsets of {Ji,..., Jgj multiplied with cumulants depending 
on subsets of {Jg,..., Jie} and B runs through the list of all blocks of the partition a. 
Now, we can use Lemma 6.1 aird insert (6.1) in (6.9). Again the largest contribution 
comes from products of eight cumulants of second order and due to = n 

if = 0 and ^ otherwise, all these combinations that occur in 

(6.9) contain a factor (5|fc.|,|fc | for i ^ j. This is either excluded in the summation in (6.8) 
or results in at least one collapsing sum, which eventually implies Var(r*^) = o(l). □ 


Proof of Remark 3.1. We only treat the variance. Under E[q, the second summand in 
(3.4) does not depend on ji,..., j 4 , so that this case vanishes thanks to 

E/ 1 (1 ~ ^iil3)(l ~ ^j2j4) \ = 

il J2,i3,l4 = l ^ ^ 

which yields the desired assertion. □ 


Proof of Theorem 3.1. Note hrst that by Theorem 6.2 and the v^nh-consistency of the 
spectral density estimators, see the proof of Theorem 2.2, we have = op(/i^/^). 

Hence, it is sufficient to prove a central limit theorem for T* — 

By using Lemma 6.2(i), we have 


ry-i* _ 


on 

■ / Kh{ui-ujk3)Kh{uj-uJk2)<^^ 

P-1 9 / 9 \ 

^ E/ E/ ( {r)^h,'^ko{r) j ^jl,mi,m2-,kilj2,mi,m2-,k2 

mi,7712=0 ji,42 = 1 \ 7-=l / 


■ ^ ' ^ra,fci ,fe2 ’ ^^2 ) 

fci ,fe2 
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with an obvious notation for w* (tJ"/!;!,7 !‘fc2) =■ fcj('^fcij^fe 2 )■ Because of 


E* 


(r) .tl (r) J 2 

\r—l 



-.-1 


|fcl| 1^21, 

i^li = 1^21, 


( 6 . 10 ) 


it follows that 

^ ^ ,^2 ^ ^ ,^2 ; ^^2 ) 

|fcl|/|^‘2] |fcl| — 1 ^ 2 ] 

=-W*+ ^ wl^^k^{-Kk^,TTk2) = W* + h-^/^p*^ + 0p,{l), 

|fcl| = |fc2] 


see Theorem 6.1 above. Setting := ■*i'fei,fc 2 ('^fei’'^fe 2 ) +'^fe 2 .fci(^fc 2 >^fei)^ we obtain 

up to a negligible term 


W* = 


E 

0<fci<fe2< L^/2J 




■IT. 


— ki,k2 


-w. 


ki , — k2 


■ w. 


— ki, — k2 


E 

0<ki<k2< [n/ 2 \ 


Since the random variables are clean by (6.10) in the sense of Definition 2.1 of De 

Jong [7], that is, £’*(1T^^ holds a.s., we are in the situation to apply Propo¬ 

sition 3.2. of his paper. Hence, for obtaining convergence in distribution (conditioned on 
the data) 

lw:^N{0,l) ( 6 . 11 ) 

El 

in probability, it remains to prove that 


Gi 

■■= E 


(6.12) 


0 <fci <fc 2 < [n/2j 



G 2 

- E 


(6.13) 


VI 

CO 

V 

V 

VI 

0 

^/2\ 



and 


G4 := E 

0 <A:i<A:2<fc3<fc4< [n/2j 

(6.14) 


are all of lower order (in probability) than r*^, that is, are all negligible. We prove the 
requested result only for the most contributing case (6.14). Because of the inherited sym¬ 
metries and the constant number of periodogram factors with frequencies 
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and uJki, it suffices to consider the representative 

^ 4 - 'y ^ (^fci ,^2 ’ ^^‘2 )^fcl ,^3 (^^1 ’ ^^‘3 )^ fc 4 , fc 2 (^^4 1 ^^2 )^ fc 4 , fc 3 (^^4 1 ^^3 )) 

0<fcl<fe2<fe3<fe4< [’’■/2J 


E 


^ Y. E r Kh{Xi-WkJK,,{X,-u;k,)dXi 


n 


0<fcl<fe2<fc3<fc4< L"/2J J8 = l ‘ 


X / Kh,{X2 — iOki)Kh{X2 — iOk;i)dX2 
J —rc 

/ 7T 

KhiXs — ujki)Kh{X^ — Wfej.) dAs 
-rc 

/ 7T 

Kh{Xi — ujki)Kh{Xi — Wfeg) dA4 
-rc 


X 


ii ,mi,m2;fei dj^ ,>Ttl ,rn2\k2^j3 ,m3,m4;fei 


X 


ji ,m3 ,m4;k3dj^ ^rn5 ,me-,ki-i ,m3 ,me-,k2 


I ^ ^ (ri) (’'l).i2 ~ 9 )■■■ 


^ri,...,7-4=1 


^ ('^•n'fc4(7-4),j7^7rfc3(r4),48 <7 ) )• 


Since the last conditional expectation above is bounded by 0 from below and by from 
above, taking the expectation above and applying (6.1) for fc = 16 results in 


E(G4) = 


E 




p-l q 

E E 


Kh{Xi — Wfci)if/t(Ai — Wfcj) dAi 


'^0<fcl<fc2<fe3<fe4<L"-/2j 7711,...,7718 = 0^1,...,48 = 1 ‘ 


X / Kh{X2 — LOk^)Kh{X2 — LOk3)dX2 
J —rc 


X / KhiXs — uJk 4 )Kh (A3 — Wfcj) dAs 
J — rc 


X J Kh{X4 —uJk4)Kh{X4 —uJk3)dX4 


xO(l). 
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Approximating all Riemann sums by their corresponding integrals and by using standard 
substitutions, the expression above becomes an 0{h) term. Similar arguments and using 
(6.1) for fc = 32 yield Var(G 4 ) = 0{hf). This completes the proof. □ 


Proof of Corollary 3.1. By comparing the results in Theorem 2.1 and Theorem 3.1, 
we have asymptotic exactness of if and only if Tq = Tq ^ holds under ilo ■ Rear¬ 

rangements of the summations over ji,j 2 ,j 3 ,j 4 in the integrands of (2.5) and (3.7) yield 
to the condition that 


(g-l)|tr(Fn(u;)Fn(u;)^)|^-- ^ |tr(Fn(u;)F,,,,(o.)^) 

^ i2,j4 = l 

1 _j, 

J2,l3 04 = 1 
il/l3 J2#i4 


(6.15) 


and 


(9-l)|tr(Fn(cu)Fn(u;)^)|^-- ^ |tr(Fn(o.)F,,,,(a.)^)|^ 

^ h,H=^ 

1 2 _T 

+ 2iq _ |i'^(^iii3(^)^i2l4('^) )l i 

ll J2,l3,i4 = l ^ ' 

3l¥=h J2#l4 


(6.16) 


have to be equal. Equalizing both quantities and another rearrangement of both last 
sums in (6.15) and in (6.16) gives the desired result. By Lemma 1 in Janssen and Pauls 
[14], this even shows the asymptotic equivalence of the tests in (3.9). Since we have 
fi — p,^ = op{h^^^) under Hq, the same holds true for the other version (/?*. 

Consistency for (p* follows as in the proof of Theorem 2.1 since T„ — con¬ 

verges in probability to -l-oo under Hi but the critical value still converges in probability 
to a quantile of a normal distribution. For the other test, we can rewrite asymp¬ 
totically equivalent as l(cj(a),oo)(2n ~ consistency follows as above since 

Tn — also converges in probability to -boo under Hi. □ 


Proof of Corollary 3.2. For q = 2 the ratios occurring in the first and in the third sum 
in (3.10) are zero and the second sum does not occur at all. In the cases (ii) and (iii), we 

_ X 

can treat | tr(Fj 4 j 3 (w)Fj 2 j 4 (a;) )p as a constant, which causes the three sums in (3.10) 
to cumulate to zero. For (iv) the right-hand side of Corollary 3.1(ii) becomes 

^(|/i2(a;)|" + |/i3(a;)|V|/23(^)|") 
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- ^(|/i2M|+ I/12MIV 23 H|^ + I/13MIV 23 M|^) 

= ^{(I/i 2M|^ - |/i3Ml^f + (I/i2H|^ - I/23MI^)^ + (|/i3MI^ - |/23M|^)^}, 

which gives the desired result. □ 

Proof of Corollary 3.3. Since r*^ converges in P ® P probability to r*^, the critical 
value c* always converges in probability to ui-q. Hence, we can close the proof 

as in Theorem 2.1. □ 

Proof of Remark 3.2. Note that the assumptions guarantee that the conditional CLT, 
see Theorem 3.1, holds true as above. Moreover, the test statistic also satisfies the CLT 
stated in Theorem 2.1 under the null, see Equations (2.4)-(2.6) in Dette and Paparoditis 
[9], which shows the asymptotic exactness under the null. Finally, consistency of the tests 
follows from Theorem 2 of their paper. □ 

Proof of Corollary 3.4. Checking through the proof of Theorem 3.1 (where the central 
limit theorem (3.3) holds under the null as well as for fixed alternatives) we see that all 
results even remain valid for the processes X”,... ,X" with spectral density given by f", 
that is, we also have convergence in distribution given the data 

in probability in this case. Since the test statistic is asymptotically normally distributed 
with mean v and variance Tq , see the proof of Theorem 2.3, the result (i) follows as in 
the proof of Corollary 3.3. 

For proving (ii), recall that our mixing conditions imply the v^nTi-consistency of the 
kernel density estimators for the processes. This shows that — /i* = op{h^^^) even holds 
under the given local alternatives (2.7) with q;„ = . Since we also have t — t* = 

op(/i^/^), the result follows from Slutsky’s lemma as in the proof of Corollary 3.1. □ 
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Supplement to “Testing equality of spectral densities using randomization 
techniques” (DOI: 10.3150/13-BEJ584SUPP; .pdf). In the supplement to the current 
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paper (cf. Jentsch and Pauly [17]), we provide additional supporting simulations for the 
asymptotic test and all three randomization tests under consideration in a variety of 
examples. 
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